Phase diagram for the z/ = quantum Hall state in monolayer graphene 



Maxim Kharitonov 

Center for Materials Theory, Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854, USA 

(Dated: January 18, 2013) 

The u = quantum Hall state in a defect-free graphene sample is studied within the framework 
of quantum Hall ferromagnetism. We perform a systematic analysis of the "isospin" anisotropies, 
which arise from the valley and sublattice asymmetric short-range electron-electron (e-e) and 
electron-phonon (e-ph) interactions. The phase diagram, obtained in the presence of generic isospin 
anisotropy and the Zeeman effect, consists of four phases characterized by the following orders: spin- 
polarized ferromagnetic, canted antiferromagnetic, charge density wave, and Kekule distortion. We 
take into account the Landau level mixing effects and show that they result in the key renormaliza- 
tions of parameters. First, the absolute values of the anisotropy energies become greatly enhanced 
and can significantly exceed the Zeeman energy. Second, the signs of the anisotropy energies due 
to e-e interactions can change upon renormalization. A crucial consequence of the latter is that 
the short-range e-e interactions alone could favor any state on the phase diagram, depending on 
the details of interactions at the lattice scale. On the other hand, the leading e-ph interactions 
always favor the Kekule distortion order. The possibility of inducing phase transitions by tilting the 
magnetic field is discussed. 

PACS numbers: 73.43.-f, 71.10.-w, 71.10.Pm 



I. INTRODUCTION 



sequence of the orbital Landau levelP^l (LLs), Fig. [l] 



Quantum Hall effects in graphen^2t2! were observed 
shortly aftei'' ^ the discovery of the material and have at- 
tracted a lot of attention ever since (for general reviews 
on graphene, see Refs. [SHHl)- Besides the "anomalous" 
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FIG. 1: Landau levels (LLs) [Eq. (l.ll] in monolayer 
graphene. Neglecting the Zeeman effect, the orbitals (i.e., 
eigenstates of the orbital part of the single-particle Hamilto- 
nian) of each LL are four-fold degenerate due to two projec- 
tions of spin and two valleys (splitting shown for illustration 
purpose and not implied). At v = filling factor, each or- 
bital of the n = LL is occupied on average by two electrons, 
while LLs with n < (n > 0) are filled (empty). The many- 
body states with identical occupation of the KK' ® s isospin- 
spin subspaces of all orbitals (only one possible occupation is 
shown) exactly minimize the energy of the Coulomb interac- 
tions, SU(4)-symmetric in the KK' (g) s space. Such states 
form the family of degenerate ground states of the i/ — 
quantum Hall ferromagnet. 



(1.1) 



characteristic of the Dirac nature of electron spectrum in 
graphene, an additional to spin, two-fold valley degener- 
acy brings in extra richness to physical phenomena [In 
Eq. (l.ll, n takes integer values, v is the Dirac velocity, 



Ib is the magnetic length dependent on the perpendic- 
ular component B± of the magnetic field, and we put 
h — 1 throughout the paper]. Among all, the n — LL 
really stands out: it is located exactly at the Dirac point 
eo = and, in each valley, K or K' , its wave-functions 
reside solely on one of the sublattices, A or B, Fig. [2] For 
n = LL, the valley KK' "isospin" is, therefore, equiv- 
alent to the sublattice AB "pseudospin" , and is referred 
to as just the "isospin" below. 

As the quality of graphene devices progressed, 
the strongly correlated quantum Hall physics clearly 
emergecP^HUl at integer and fractional filling factors v. 
To date, one of the most intriguing questions concerns 
the nature of the v — quantum Hall state, in which 
the orbitals of the n = LL, four-fold degenerate in the 
KK' (g) s isospin-spin space in the absence of the Zeeman 
effect, are occupied on average by two electrons. Fig. [T] 
The interest is largely motivated by the strongly insu- 
lating behavior o f the state, initially observed in sa mples 
on Si02 substrat d^^ l ^"^ ! later - in suspended sampleJ^^E^, 
and quite recently - in samples on boron nitride sub- 
strate^. The tendency of the higher quality samples to 
be more resistive signifies that the insulating behavior is 
an intrinsic property of an ideal defect-free system, rather 
than it is due to disorder-induced localization effects. Al- 
though the interaction-driven character of the insulating 
v = state is apparent, its precise nature remains an 



2 




FIG. 2: Graphene honeycomb lattice and the structure of 
the wave-functions of the n = LL at the lattice scale. In 
each valley, K or K' , the wave-functions reside on only one 
sublattice, A (left) or B (right). 

open challenge in the graphene field. 

A l ot of th eoretical activity has been Ae- 
^o^g jMMIiaflllllMl to the properties of the ly = 
quantum Hall state in graphene (see also Refs. 139)401 
for reviews). Due to the quenched kinetic energy, in 
sufficiently clean samples and/or high enough magnetic 
fields, electrons in a partially filled LL present a strongly 
interacting system. While at arbitrary (e.g., fractional) 
filling factors u such systems pose a formidable theoreti- 
cal challenge, an appealing property of the integer fillings 
is that, in the leading approximation, the family of the 
many-body ground states can be found exactly. For all 
integer v in graphene, including v = Q, this solution 
is provided b y the gen eral theory of quantum Hall 
ferromagnetisrrpSElEllHsl (QHFM) for multicomponent 
systems, in which spin and valley degrees of freedom are 
united into one SU(4) "spin" . (Alternative views on the 

V = Q state based on the idea of "magnetic catalysis" 
can be found in Refs. I34H381 ) 

The idea of constructing the exact solution can be 
traced back to the Hund's rule in atomic physics: specif- 
ically for the V = Q state, the energy of the repul- 
sive Coulomb interactions is minimized by the many- 
body states, in which the four-dimensional (4D) KK' ®s 
isospin-spin subspace of each orbital of the n = Q LL is 
occupied by two electrons in exactly the same way, Fig.[T] 
Two key provisos that make this result exact are (i) ab- 
sence of the single-particle or many-body perturbations 
that break the SU(4) symmetry in the KK' ® s space 
and (ii) absence of the interaction-induced electron tran- 
sitions between different LLs (LL "mixing"). For SU(4)- 
symmetric Coulomb interactions, the ground state is de- 
generate with respect to the choice of occupation of the 
isospin-spin space and there is no preference between or- 
dering of the spin vs. isospin degrees of freedom. 

This exact result is the cornerstone of the theory of the 

V = Q state. In a real graphene system, neither of the 
conditions (i) or (ii) is satisfied precisely, and the result 
is generally not exact. However, it is possible to develop 
the low-energy quantum field theory for the v = Q quan- 
tum Hall ferromagnet (QHFM), in which the deviations 
from these conditions are consistently taken into account. 



As we show here, for Coulomb interactions of moderate 
strength /v ^ 1, as is the case in graphene, the ef- 
fects of LL mixing can be systematically taken into ac- 
count in the large- iV approximation. Further, the factors 
that break the symmetry in the AB (g) s space typically 
have smaller energy scales than the Coulomb interactions 
and can be taken into account perturbatively. They are, 
nonetheless, extremely important, as they lift the rich 
degeneracy of the ground state favoring certain orders, 
the physical properties of which may differ substantially. 

The simplest single-particle mechanism is the Zee- 
man effect, which breaks the spin symmetry and nat- 
urall y favor s the fully spin-polarized ferromagnetic (F) 
stat e! 1 ^'^ *^. Since spin-orbit interaction is weak in 
graphene, the Zeeman effect is practically the only rele- 
vant factor that affects the spin symmetry, while other 
key perturbations break the isospin symmetry. An anal- 
ogous regular "Zeeman" field for the KK' isospin is vir- 
tually impossible to implement in a controlled way in 
real graphene, although a random one can arise from the 
short-range disorder. 

However, even in an ideal defect-free sample electron- 
electron (e-e) and electron-phonon (e-ph) interactions 
necessarily break the valley and sublattice symmetry 
at the lattice scales^^ Having many-body ori- 

gin, these mechanisms give rise to the isospin aniso tropy 
rather than the Zeeman- type fields. Existing studie^^HHl 
of the lattice effects of e-e interactions on the v = Q state 
in the framework of the QHFM theory were carried out 
using the tight-binding extended Hubbard model, with 
adjustable interactions at the lattice scale and Coulomb 
asymptotic at larger distances. Depending on the in- 
teractions at the lattice scale, the competition between 
the F and charge-density-wave (CDW) ground states was 
predicted in Ref . 1221 while the numerical mean-field anal- 
ysis of Ref. [29] predicted cither the CDW or the antifer- 
romagnetic (AF) ground state. Electron-phonon inter- 
actions, on the other hand, were predicted to favor the 
fully isospin-polarized states with either the Kekule dis- 
tortiorpMU (kd) or CDW^I orders. We also mention 
that similar phases were predicted in Ref. [37] within a 
different framework of "magnetic catalysis" . 

The variety of the proposed ground states poses a ques- 
tion whether the above list is exhaustive, i.e., whether, 
for a given source of the isospin anisotropy, one can find 
all possible orders that could be realized in a general case 
scenario. Attempting to answer this question, in this 
paper, we perform a systematic analysis of the isospin 
anisotropics of the v = Q QHFM arising from the short- 
range e-e and e-ph interactions, without appealing to 
any specific lattice model. Starting from the most gen- 
eral form of e-e interactions allowed by symmetry in the 
Dirac Hamiltonian and taking the leading e-ph interac- 
tions into account, we derive the low-energy QHFM the- 
ory in the presence of the isospin anisotropy and obtain 
a phase diagram for the v = Q state. The diagram ob- 
tained by separately minimizing the energy of the generic 
isospin anisotropy (i.e., considered without any assump- 
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tions about the nature or properties of the underlying 
interactions, whether e-e or e-ph, repulsive or attrac- 
tive) consists of four phases: SU(2)-spin-degenerate F 
and AF phases, and spin-singlet CDW and KD phases 
with Z2 and U{1) isospin degeneracies, respectively. In- 
cluding the Zeeman effect does not alter the CDW or KD 
phases, but removes the degeneracy of the F phase and 
transforms the AF phase into a canted antiferromagnetic 
(CAF) phase with noncoUinear spin polarizations of sub- 
lattices. 



We also consider the critical renormalizations of the 
isospin anisotropy by the long-range Coulomb interac- 
tions and address the question whether one can rule out 
certain states from the phase diagram based on the re- 
pulsive nature of the Coulomb interactions and attrac- 
tive nature of the phonon-mediated interactions. We 
arrive at an important conclusion that the short-range 
e-e interactions could favor essentially any state on the 
generic phase diagram: F, CAF, CDW, or KD. The rea- 
son for this are peculiar properties of the renormaliza- 
tions, which allow for sign changes of the e-e coupling 
constants, switching the interactions from repulsive to 
attractive in certain valley-sublattice channels. As a re- 
sult, unless a reliable numerical estimate for the bare 
short-range e-e coupling constants is provided, one can- 
not theoretically rule out any possibility. In contrast, 
the phonon-mediated interactions remain attractive un- 
der renormalizations and the leading e-ph interactions 
always favor the KD order. 

Among potential practical applications of the present 
work, we also study the transitions that could be induced 
between the obtained phases by the tilting magnetic field. 
Once combined with a more details analysis of the charge 
excitations of these phases, to be presented elsewhere^'^, 
our findings could be used to identify the particula r phase 
realized in the experimentally observed ^ 1 5 1 1 7 | i s] insulat- 
ing ly = state. 

The rest of the paper is organized as follows. In Sec.|TTj 
the low-energy Hamiltonian, which describes electron dy- 
namics in the vicinity of the Dirac point, is presented, 
the basic properties of Landau levels in graphene are dis- 
cussed, and the projected Hamiltonian for the n — 



Landau level is derived. In Sec. Ill the low-energy quan- 
tum field theory for the — quantum Hall ferromagnet, 
which includes the effects of the isospin anisotropy, is de- 
rived. In Sec. |IV[ the Landau level mixing effects are 
considered in the large-TV approximation and the criti- 
cal renormalizations of the isospin anisotropics are stud- 
ied. In Sec. |V] the phase diagram of the ly — quantum 
Hall ferromagnet in the presence of the isospin anisotropy 
and Zeeman effect is obtained. The possibility of in- 
ducing phase transitions by tilting the magnetic field is 
discussed. Concluding remarks and connection to the 



II. MODEL AND HAMILTONIAN 

We start the analysis by writing down the low-energy 
Hamiltonian 



H — Ho + He- 



c-phj 



(2.1) 



which describes electron dynamics in graphene in the 
vicinity of the Dirac point. The terms Hq, -ffc-ci £^nd 
Hc-ph, describing noninteracting electrons, e-e and e-ph 
interactions, respectively, are discussed in the next three 
subsections. Our Hamiltonian and the choice of basis are 
identical to those of Refs. 133133 with some differences in 
notation. 



A. Basis and single-particle Hamiltonian 

At atomic scales, the single-particle electron Hamilto- 
nian can be written as 



-Ho = / d'r^Ur) 



2m 



U{r) 



*.(rO. (2.2) 



Here, is the electron field operator, f— {x,y,z) is 

a continuous three-dimensional (3D) radius-vector, and 
a i is the spin projection (summation over a is im- 
plied). The self-consistent periodic potential U{r) of the 
graphene honeycomb lattice constrains electrons around 
the 2 = plane and has a Cg^, point group symmetry 
within the plane. This symmetry dictates the follow- 
ing properties of the graphene band structure. Exactly 
at the Dirac point, taken to be at zero energy e = 0, 
there are four orthogonal Bloch-wave solutions uxa{^j 
UKB{r), UK'Ai'r), UK'sif) of the Schrodinger equation 
associated with Eq. ( |2.2[ ). The indices K and K' re- 
fer to different valleys in the Brillouin with wave vectors 
K = ;^(ai— a2) and K' — — K, respectively (ai 2 are the 
primitive translations of the honeycomb lattice, shown in 
Fig. and oq = |ai_2| ~ 2.46A is the lattice constant), 
and the indices A and B indicate that the wave-functions 
are predominantly localized at the positions of the A and 
B sites. The solutions corresponding to different valleys 



are related as UK'A{r) — "^^^(r) and UK'sir) 



'■KB 



(r). 



experiment are presented in Sec. VI 



For the excitation energies e much smaller than the 
bandwidth, one may expand the electron field in terms 
of the e = solutions, 

+ i'K'Acr{r)uK'A{r) + 1pK'Ba{r)uK'B{r}{2-3) 

The Dirac field operators V'Ao-(r), A = KA, KB, K'A, 
K'B, are functions of a two-dimensional (2D) continu- 
ous radius vector r = (x, y) and vary at scales much 
larger than the atomic one ag. To ensure the standard 
normalization of the fields reflected in anticommutation 
relation 

{^xa (r) , Vl- V (r')} = 5xy <5..' S{r - r') , 
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where 5{r — r') is a 2D delta function at large scales 
and {, } is the anticommutator, the Bloch wave- functions 
ux{r) must be normalized as 



d^ful{r}uy{r) - 3|[ai x asH^AV, (2.4) 



3uc 



where the integration is performed over the tripled unit 
cell (3uc), which contains six atoms. 

The Dirac fields can be joined in a vector as 



(2.5) 



KK'(g)AB 



The advantage of this ordering is that it gives the 
most symmetric representation of the Dirac Hamiltonian. 
Since this way the sublattice indices in the K' valley are 
interchanged, to avoid c onfu sion, we denote the sublat- 
tice space of the basis (2.5) as AB, and the whole 4D 
space ~ as the direct product KK' (g) AB. With spin 
included, the low-energy electron degrees of freedom are 
described by the eight-component field operator 



^(r) 



(2.6) 



in the direct product KK' <Si AB <Si s of the valley {KK'), 
sublattice (AB), and spin (s) spaces. 

The symmetry properties of the Bloch wave-functions 
u\ (r) at the Dirac point are sufRcient to derive the many- 
body low-energy Hamiltonian in the basis of The 
single-particle Hamiltonian, obtained from Eq. (2.2), has 
the form 

Ho^ U^vij\v) V J2 Toc.{pa--^A^)-ezS, ^j{r), 

L a—x,y 

(2.7) 

where Pa — — iV^, V — {dx,dy), and v ~ lO^cm/s is 
the velocity of the Dirac spectrum. Here and below, for 
a,/3 = 0,x,y,z, 

with the unity (tq — 1) and Pauli {tx, Ty, Tz) matrices in 
the corresponding 2D subspaces. In Eq. ( |2.7[ ), we intro- 
duced the orbital and spin effects of the magnetic field 
[not written in Eq. (2.2|], described by the vector poten- 
tial Aa{r), rot A = (0, 0, B±), and the Zeeman term with 

B^, and 



ez = ^J.BB, B = 



Bl 



Sz 



1 



KK' 



1 



AB 



We assume arbitrary orientation of the total magnetic 
field B — (By, 0, i?_L) relative to the plane z = of 
graphene sample; the z direction in the spin space points 
along B and is not necessarily perpendicular to the sam- 
ple. 



B. Electron-electron interactions 

The most general form of the spin-symmetric e-e inter- 
actions in the low-energy Hamiltonian (2.1 ) can be writ- 
ten dowiP^ solely based on the symmetry considerations 
as 



'o-o.O 



(2., 



Here, 



^o-c,o = 2 y d^i-d^r' [V'^(r)^(r)]^o(r - r')[V;^(r')^(r')] 

(2.9) 

describes the long-range Coulomb interactions, Vo(r) — 
e^/|r|, symmetric in valley-sublattice space KK' (E) AB 
and 
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d\J2 9c.p[^^ir)%.f!i'ir)? (2.10) 

a.,l3 



describes the short-range e-e interactions that break the 
valley and/or sublattice symmetry. The summation 

includes all combinations a, /3 = 0, x, y, z, of the 

valley a and sublattice /? channels, except for the sym- 



metric one a — /3 — 0, which is given by Eq. (2.9 1. In 



Eqs. (2.9) and (2.10), and below, normal ordering of the 



operators is understood. 

The symmetry of the honeycomb lattice yields the fol- 
lowing relations between the coupling^^, 

.9-L-L — 9xx 9xy — 9yx — 9yy^ 



'Jz± 



Qzyi 



= .92:0 — 9yVn go± = dOx — QOy 

Thus, the asymmetry of the interactions in KK' (g) AB 
space is described by eight independent coupling con- 
stants g_L_L, g_i_z, gz±, g^z, g±o, 9z0, 9o±, goz- 

Although, of course, the origin of both symmet- 
ric [Eq. (2.9)] and asymmetric [Eq. (2.10)] e-e interac- 
tions are the actual Coulomb interactions, for brevity, 
we will refer to them as the "Coulomb" and "short- 
range/asymmetric e-e" interactions, respectively. 

To lowest orders, the expressions for the coupling con- 
stants g^p can be obtained by considering the Coulomb 
interactions in the atomic-scale model. 
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dVdV^1,(f)^'T,(r' 



*.'(r=')*.(r). 

(2.11) 



Substituting the expansion (2.3) for ^cr(r) into Eq. (2.11 ) 
and using the slow variation of V'(r) at atomic scales, one 
obtains Eqs__(p^, ( pj] ), and ( |2.10[ ) with the first-order 
expressions^ 



9al3 



dV 



3uc 



3|[ai X as] 



I Pa/3 (t^) 

(2.12) 



for the couphngs gap- In Eq. (2.12) 



are the real densities in a given vallcy-sublattice channel 
aB and 



UKA{r) 
UKB {r) 
UK' sir) 

-Uk'a{t^ , KK'<^AB 



From the orthogonality properties (2.4) of the Bloch 
wave-functions, it is clear that the integrands in 
Eq. (2.12), as functions of r — f*, decay over several unit 



cells and the asymmetric interactions are indeed short- 
ranged. This emphasizes the fact that breaking of the 
valley-sublattice symmetry arises from atomic scales. 

Using the relations u*j^iji^(r) = uxAif) and u*j^,g{r) = 
UKsi'T), we see that the densities /9oa(r) — Paa{r) = 0, 
a = x,y,z, vanish identically. Therefore, g^^^ — g^^J = 

9±o = gio = 0' a-nd the couplings ga±, g^z, g±o, 9z0, 
although not prohibited by symmetry, vanish in the first 

order. The nonvanishing expressions g^^ for gap-, with 
a = or /3 = 0, arise in the second order in the 
Coulomb interactions and involve virtual transitions to 
other bands; we do not present these expressions here. 

For future discussion in Sees. IV and [V] we note the 
following properties. The first-order microscopic expres- 
sions (2.12) for the coupling constants have the form of 



the electrostatic Coulomb energy for the density distri- 
butions Pap (r) - Since it is well knowrP^l that the elec- 
trostatic energy is positive- definite, i.e., positive for any 
nonvanishing charge distribution, we conclude that all 
nonvanishing first-order expressions (2.12) must be pos- 
itive, 

9±{ > 0, g^2l > 0, g'-'l > 0, g^^J > 0. (2.13) 

On the other hand, the second-order expressions have to 
be negative. 



J2) 



< 0. 



(2.14) 



One can expect the lowest-order expressions 



(2) 

[Eq. (2.12)] and to provide accurate estimates for the 
couplings gap in the limit of weak Coulomb interactions, 
e^/w <^ 1- For stronger interactions, e^/v ~ 1, as in real 
graphene, this is not necessarily the case. The reason is 
that the short-range interactions renormalize themselves 
at energies on the order of bandwidth w/ap. As an illus- 
tration of this fact, the diagrams of the low-energy Dirac 
theory involving just the short-range interactions contain 
ultraviolet divergencies: schematically, each extra order 
produces a relative factor gap J^^°'° dq/v ~ gap/ivao) ^ 
/v- These renormalizations change the magnitude and, 
possibly, the signs of the couplings gap in certain chan- 
nels. 



gap 




FIG. 3: In-plane optical phonon modes with the strongest e- 
ph coupling, (left) Linear combination of the two degenerate 
E2 modes {'u,x{r),Uy{r)) with the phonon wavevector at F 
point, (right) Linear combination of degenerate Ai, B\ modes 
{ua{v),Ub{Y)) with the wavevector at K,K' points. 



Thus, it is more reasonable to treat the couplings gap 



in Eq. (2.10) as the bare inputs of the low-energy the- 
ory, without any assumptions about their signs and rela- 
tive values, and consider all possibilities. This is the ap- 
proach we choose in the of paper. An order-of-magnitude 
estimate for the bare couplings, valid for both weak 
{e^/v ^ 1) and moderate (e^/w ~ 1) Coulomb inter- 
actions, is 



9aP 



(2.15) 



Anticipating the results of the next sections, the 



asymmetric short-range e-e interactions (2.10), although 
weaker than the symmetric Coulomb ones (2.9), appear 
to be play a crucial role in the physics of the v = 
QHFM. Possible relations between and the signs of the 
couplings gap become especially important, as they de- 
termine the properties of the isospin anisotropy and, as 
a result, the favored ground state order. The implica- 
tions of the potential sign restrictions on gap, suggested 
by Eqs. 
andM 



C. Electron-phonon interactions 



(2.13) and (2.14), wiU be discussed in Sec. IV 



Besides the short-range e-e interactions (2.10), another 
source of the isospin anisotropy in the 1^ = QHFM 
comes from e-ph interactions. Electrons in graphene 
couple most efficiently to the following in-plane optical 
phonons: two E2 modes with the phonon wave- vector at 
the F point and Ai,Bi modes with wave- vector at K, K' 
point (following the classification of Ref. shown in 
Fig. [3j The corresponding e-ph interactions can be de- 
scribed by the Hamiltonian 



H, 



c-ph 



d^rVj^(r){F£;2 [TzyU.^{r) - TzxUy{r)] 



FAiiTxzUair) +TyzUb{r)]}^{r)- 



(2.16) 



The two degenerate E2 modes have the frequency uje2 ~ 
0.196eV. The Ai and Bi modes are also related by sym- 
metry and have the same frequency ujai ~ 0.170eV and 
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same coupling constant Fa^ = -Fbi • The lattice deforma- 
tion due to Ai , Bi modes is that of the Kekule distortion. 
The two-fold degeneracy of the modes allows for arbitrary 
in-plane displacement of a given atom in Fig.js] (left) and 
(right), while the displacements of the remaining atoms 
in a tripled unit cell are related by the corresponding 
symmetry. 

The phonon dynamics is described by the correlation 
fimctions of the displacement operators ■u^(r), 



So 



2Mu, 



-Sir) 



2^ 



(2.17) 



in the Matsubara representation. Here, M is the mass of 



the carbon atom and sq 



Cq is the area per carbon 



atom. We will perform calculations at zero temperature 
only, in which case a; is a continuous frequency. 
The order-of-magnitude estimate 



(2.18) 



for the coupling constants follows from the dimensional 
analysis of Eq. (2.161 and the electrostatic origin of e-ph 
interactions. 



D. Landau levels in graphene 

In this section, we briefly present the single-particle 
basis of the problem and emphasize the key properties of 
the n = LL. 



Solving the Dirac equation associated with Eq. (2.7) 
in the gauge A(r) = (0, i?_LX, 0), in the AB sublattice 
space of each valley K and K\ one obtainJ^^li' the wave- 
functions 



for all integer n 7^ and 



AB 



(r|Op) 



_ ( 4>q{x - Xp) 





AB 



(2.19) 



(2.20) 



with the orbital energies e„ given by Eq. (1.1 1. Here, 
4>\n\{x) are the harmonic oscillator wavef unctions, Xp ~ 
pPg is the "guiding center", and Ly is the size of the 
sample in the y direction, introduced to discretize the 
momentum quantum number p. 

The complete set of the single-particle eigenstates in 
the KK' (g) AB ® s space is given by 



with 11 = K, K' , G =t,i, and 



(2.21) 



\K) = 



KK' 



lt> = 



,u> = 



The fi eld operator (2.6 1 can be expanded in the ba- 
sis (|2.21|) as 



■0(1") = X! ^"('')' '^«('') ^^A"n'VU-o)Cnvyi.a. (2.22) 



where Cnp^a are the annihilation operators. 

The n — LL with eg = is located exactly at the 
Dirac point and possesses arguably the most peculiar 
properties: in each valley, K or K', its wave-functions 



(2.21 ) resides solely on one (actual) sublattice, A or B, re- 
spectively [see Fig. [2j Eq. (2.20), and recall the accepted 
ordering (2.6)]. Therefore, for each spin projection a, the 
part 



'0o<T(r) = ^{ra\Opiia) 



( iJ'OKAair) \ 



ffjQK'Bair) 





/ KK'®AB 

(2.23) 

of the field operator (2.22) pertaining to n = LL 



has only two nonvanishing components, li^oKAai^) and 
VW'Bcr(r), whereas 



^OKB<T(r) = V'OK'ylcr(r) = 0. 



(2.24) 



Thus, for n = LL the valley [KK') and sublattice {AB) 
degrees of freedom are essentially equivalent, K ^ A, 
K' ^ B. Further on, when discussing n = LL below, 
we refer to this 2D degree of freedom {ipQKAa-,'4'0K'Ba-) 
as just the ''KK' valley isospin". Accordingly, we join 
the nonvanishing components of the n = LL operator 
(|2.23|) in a 4D vector 



V'o(i-) 



f ■0oifAt(r) ^ 
■foKAiir) 

■0oK'Br(i") 

V ^0K'B;(r) / 



(2.25) 



KK't 



in the isospin-spin space KK' (g) s. 



E. Projected Hamiltonian for n = LL 

When addressing the many-body aspects of the i/ = 
state, as a starting point, one may neglect the contribu- 
tions from n ^ LLs and restrict oneself to the dynamics 
within n = L L, de scribed in terms of the field (2.25). 
From the form (2.23) of ipo{r), we obtain 



_ i i'lTalpO, ^ = 0,2, 

0, P^x.y, 



for a = 0,x, y, z, where 



KK' 



KK' 



0,x,y, z. 



(2.26) 



(2.27) 
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are the KK' -isosp'm matrices. I.e., electrons in the n = 
LL couple directly only to the source fields with Tao 
and Taz vertex structure in KK' (g) AB space. This is 
a consequence of the properties of the n = LL wave- 
functions. 



Sub stituti ng -0 (r) in the form ( |2.22[ ) into Eqs. ( |2.9| , 
(2.10), and (2.16) and retaining only the n — LL com- 



ponent (2.23), we obtain the bare projected Hamiltonian 



in terms of the field (2.25) 



(0) 



e-e,0 



(0) 



H, 



(0) 



c-ph ' 



o-(O) 
^0 



(2.30) 

Hi'l = lld\ 9c.[i'Ur)TMr)r, (2.31) 

a—x.y.z 

with 5a = g„o + cjaz, and 



6z / d'rijlir)SzMr), Sz = 1 



(2.28) 



(2.29) 



C-ph 



7^U6(r)]V;o(r). 

(2.32) 

In the single-particle Hamiltonian since the ki- 

netic energy eq = 0, only the Zeeman t erm i s present. 



In and H^^^, due to the property (|2.26|), only the 



short-range e-e interactions with 9^0 and g^^, a — x,y, z, 
couplings and e-ph interactions with Fa^ and Fb^ cou- 
plings remain. The coupling goz, although also does not 
vanish, produces a symmetric term (x (7oz[V'o(r)'0o(r)]^, 
which may be neglected compared to the Coulomb part 
(2.30). We also mention that the trigonal warping ef- 
fect cx Txz,Tyz does not couple to n 
perturbative level. 



states at the 



III. f^O QUANTUM HALL FERROMAGNET 
IN GRAPHENE 

A. Basic concept and exact result 

At integer filling factors v, interacting multi- 
component quantum Hall systems are described by 
the general theory of the quantum Hall ferromag- 
netisrrPMIIilHini, Jtg central point is that, as long as 
electron dynamics may be effectively restricted to the 
partially filled LL (sufficient conditions for this will be 
discussed in Sec. IV) and the interactions are symmet- 
ric in the "spin" space of discrete degrees of freedom, the 
family of the many-body bulk ground states can be found 
exactly as follows. In order to minimize the energy of 
the Coulomb repulsion, one makes the orbital part of the 



wave-function totally antisymmetric, thus putting elec- 
trons, for a given density, as far apart from each other as 
possible. Since one has on average an integer number of 
electrons per orbital, such wave-function can be realized 
if electrons occupy the discrete states of all orbitals in 
exactly the same fashion. 

Specifically for the v — state in graphene, which 
hosts two electrons per four-fold degenerate orbital of 
the n = LL (Fig.[l]), the many body wave-function can 
be written as 



n E' 

p \\a,\'a' 



\(T,\' a' 
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ID). (3.1) 



Here, |0) is the "vacuum" state with completely empty 
n > Q LLs and completely filled n < LLs. Each fac- 
tor in the product J^^ creates a pair of electrons in the 
state <i> = {^\cr,\'a'} (A, A' = A, B an d a, a' =t,i) at 
orbital p of the n = LL, see Eq. (2.21 ) and we identify 



K A and K' o B; the antisymmetric two-particle 
spinor $ describes the occupation of the 4D KK' (g) s 
isospin-spin space of eac h or bital by two electrons. The 



sum 



\cr,X 

off-diagonal elements of 



in Eq. (3.1) goes over the upper-right 
$ and we normalize the spinor 



according to the number of particles per orbital. 



E 

Xa.X'a 



1$ 



Act, A' 



(3.2) 



For the purpose of illustrating the said exact property, 
in this subsection, we simply neglect the other n ^ LLs. 
To make this justified, one may temporarily assume here 
that the Coulomb interactions are weak, e^/v 1, and 



in Sec. IV we demonstrate how the LL mixing effects can 
systematically be taken into account for stronger inter- 
actions e'^ /v ~ 1 within the large- A'^ approximation. 

Acting with the Coulomb interaction Hamiltonian 
(2.30), symmetric in KK' (E) s space, on the state (3.1), 



we obtain that ^' is an eigenstate, 

if and only if $ satisfies the constraint 

(3.3) 

The energy of the state equals 



^0 = jE[4^o'(^,fc')-2Kr(^,fc') 



(3.4) 



k,k' 



where 
F^(fc,A:') = 

V%k,k') 



1 

L,, 



dqx 

2ti 



'yQ{(lx,qy = 0), 



dgx 
27r 



'Vo{qx,k-k') 
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are the "direct" (Hartree) and "exchange" (Fock) ma- 
trix elements, respectively, and Vb(q) = 27re^/|q|, q = 
{qx,%), is the Fourier transform of the Coulomb poten- 
tial. The energy Eq can easily be calculated explicitly. 
Equation (3.3 1 is a necessary and sufficient condition 



for the two-particle spinor $ to be a Slater-determinant 
state, 



XaOXb-XbO Xa, 



(3.5) 



described by two orthogonal single-particle spinors Xa.b 
in KK'®s space; the symbol o denotes the direct product 
of KK' spaces of two electrons. We see that not every 
antisymmetric spinor $ delivers a many-body eigenstate 
of the interaction Hamiltonian H^^^^q: for example, the 
spin-singlet isospin-triplet state with zero isospin projec- 
tion is not an eigenstate. However, a ny s tate (3.1 ) with <i> 



in the form of a Slater determinant (3.5) is an eigenstate 
with the energy Eq and the ground stat^^ is, therefore, 
degenerate. 

Let us count the number of degrees of freedom parame- 
terizing the ground state 4*. An arbitrary antisymmetric 
spinor $ has (six complex) = (twelve real) degrees of free- 
dom. Fixing its norm and inconsequential overall phase 
factor leaves ten real parameters, and the complex con- 
straint (3.3) reduces this to the final eight real parame- 
ters. 

Each Slater-determinant state <& is uniquely specified 
by the 2D subspace of the KK' (g) s space, occupied by 
two electrons and generated by the vectors Xa.t- This es- 
tablishes a one-to-one correspondence between the states 
(3.5) and the elements of the Grassmannian manifold 



Gr(2, 4), known as the Pliicker embedding in mathemat- 
ical literatur^^ the constraint (3.3 1 is called the Pliicker 



relation. The parametrization of the occupied subspaces, 
generated by Xa,b, and therefore of the Grassmannian 
Gr(2,4), is efficiently realized by the matrix 



Using Eqs. (3.1) and (3.51, this gives 



P = XaXl + Xbxl 



(3.6) 



(3.7) 



in the matrix form in the KK' <Si s space. The single- 
particle density matrix P satisfies the properties of a 
hermitian projection operator 



pt = p p2 ^ p 

and also, for the doubly-filled 1^ = state, 
trP = 2. 



(3.8) 



(3.9) 



The matrix P [Eq. ( 3.6 )] plays the role of the order pa- 
rameter of the broken-symmetry state (3.1); the observ- 



ables and coupling of the state (3.1 ) to various perturba- 
tions can be expressed through it. The matrix P is re- 
lated to the matrix Q {R in Ref . |4T|) . commonly usecP^'^ 
in the QHFM theory, as P = i(i + Q). 



Speaking of symmetries, since ^ [Eq. (3.1 1] is a ground 
state for any Slater-determinant state (3.5), any SU(4) 
transformation in the single-particle KK' ®s space keeps 
the energy Eq invariant. However, any SU(2) x SU(2) x 
U(l) transformation, corresponding to independent ro- 
tations within the subspace of the occupied states, gen- 
erated by Xa and Xfcj and its orthogonal complement - 
the subspace of empty states, not only does not change 
the energy E'o, but also leaves the state $ intact. There- 
fore, t he sy mmetry of the v = Q QHFM state 4" is de- 
scribecPJ^by the factor group SU(4)/[SU(2) x SU(2) x 
U(l)] = U(4)/[U(2) x U(2)]. This group also determines 
the transformation properties of the order parameter P. 
The dimensionality of the space of matrices P as an 
U(4)/[U(2) X U(2)] manifold is 4^ - 2^ - 2^ = 8, which 
agrees with the number of the physical degrees of freedom 
of the Slater-determinant states $ [Eq. (3.5)]. 



B. Energy functional of the = QHFM 

The exact result of the previous section lays down the 
basis of the QHFM theory. In the presence of the isospin- 
asymmetric interactions ( |2.31 ) and (2.32), the state (3.1) 
will generally no longer be an exact ground state. Be- 
sides, it is desirable to know not only the ground state of 
the system, but also the excitations. Provided the energy 
scales (per orbital) of these perturbations and excitations 
are small compared to the energy mi n(^, j ^^) of the 
screened Coulomb interactions (see Sec. IV A below), the 
local deviations of the actual many-body eigenstate from 
the QHFM state (3.1) are also minor. This makes pos- 



sible to develop a systematic low-energy quantum field 
theory that describes the dynamics of the system. Such 
theory has the form of a U(4)/[U(2) x U(2)] sigma-model 
for the order parameter P{t,r), which acquires time and 
coordinate dependence. 

The rigorous derivatiorPlEll of the sigma-model in- 
volves a procedure of projecting onto the QHFM state 
with a given order para meter P{t,r). Proceeding along 
the standard stepJ^^^, we arrive at the following La- 
grangian of the v = QHFM 



L[P{t,r)]^K[P{t,r)]-E[P{t,r)]. 



(3.10) 



Here, K[P{t,r)] is the kinetic term containing the time 
derivative of P(t, r); it is most simply expressed in terms 
of the single-particle spinors Xa,6(i, i") [Eq. (3.7|] 



K[P{t,r)]=i I £^{xidtXa + xldtXb 



The energetics of the v = Q QHFM is described by the 
energy functional 



E[P{t,v)] 



d\ 

27r/| 



[8o{P) + £^{P) + £z{P)]- (3.11) 



Here, 



£^[P) = p^tr[VPVP] 



(3.12) 



is the gradient term characterized by the stiffness ps and 

£ziP) = -^ztiiS,P] (3.13) 

is the Zeeman term characterized by the energy ez — 
fisB. Most importantly, 



a—x,y,z 



t^{P) = iv[Tc,P] tr[r„P] - tv[%P%P] 



(3.14) 



(3.15) 



is the isospin anisotropy energy arising from the short- 
range e-e [Eq. ( |2.8[ )] and e-ph [Eq. (2.16)] interactions, 
asymmetric in the valley-sublattice space. The KK'- 



isospin matrices Ta were introduced in Eq. (2.27). As we 



will see below, due to symmetries of the e-e and e-ph cou- 
pling constants, the anisotropy energies Uq. for a = x,y 
isospin channels are equal. 



(a) 




(b) 








FIG. 4: Diagrammatic representation of the terms in the 
energy functional (3.111 of the i/ = QHFM, defining the 
bare (lowest order in interactions) values of parameters, (a) 



Diagram for the gradient t erm £ o(P) [Eq. (3.12 1], determining 
the bare stiffness [Eq. ( 3.16 1] ; the wavy line stands for the 
Coulomb interaction [Eq. (2.30)]. (b) an d (c) Diagr ams fo r 
the isospin anisotropy term £<>{P) [Eqs. (|3.14 \ and ( 3.15 \] 



determining the bare anisotropy energies [Eqs. (3.171 



(3.181, and (3.19l]. The dashed line represents either the 
short-range e-e [Eq. ( |2.31 1] or e-ph [Eq. ( |2.32 1] interactions. 
The Hartree (b) and Fock (c) contributions produce the first 
and second terms in Eq. (3.151, respectively. 



Uj_ — Ux = Uy. 

Thus, the isospin anisotropy is fully characterized by two 
energies, u±^ and u^. 

The expressions (3.14) and (3.15) for the isospin 
anisotropy energy of the v ~ Q QHFM constitute 
one of the key results of the present work. This is 
a generic form of the anisotropy, arising from spin- 
symmetric two-particle electron interactions (including 
phonon-mediated interactions) with arbitrary structure 
in the valley-sublattice space. 

In Sec. |V) we minimize the energy functional (|3.11|) 
to obtain a phase diagram of the v = Q QHFM in the 
presence of the isospin anisotropy and Zeeman effect in 
the space of parameters (uj_,Uz,e2). In the rest of this 
and in the whole next section, we discuss the expressions 
for the stiffness Ps and anisotropy energies mj_ in terms 
of the microscopic parameters of the Dirac Hamiltonian. 
Their bare values p^s^^ and u\^^, obtained in the lowest 
order in interactions, can be determined from the pro- 
jected Hamiltonian, Eqs. ( |2.28P -( |2.32[ ). 

The gradient term £o{P) arises from the Fock free- 
energy diagram in Fig. |4][a) , in which one needs to take 
the spatial inhomogeneity of order paramete r P(t, r) into 
account. This yields the standard expressiorPilSIlIl 



Ps 



1 



16^2^ Ib 



(3.16) 



for the bare stiffness. 

The isospin anisotropy term £o{P) can be represented 
by the free-energy diagrams in F igs. [4|^b) and|4](c); the 
first and second terms in Eq. ( 3.15 ) arise from the Hartree 



(b) and Fock (c) diagrams, respectively. The diagrams 
for the short-range e-e [Eq. (2.31)] and e-ph [Eqs. (2.17) 
and (2.32)] interactions have the same form and the 
dashed line in the figures stands for either the short-range 
e-e or e-ph interactions. This way, for the bare anisotropy 



energies in terms of the valley-sublattice asymmetric cou- 
plings, we obtain 



,(c-c,0) 



,(o-ph,0) 



a =_L, z, 



where 



and 



,(o-o,0) _ 



(e-ph,0) 



(.9aO + gaz), a =_L, Z, 



J±z 



^(e-ph,0) ^ 



(3.17) 



(3.18) 



(3.19) 



are the anisotropy energies due to short-range e-e and 
e-ph interactions, respectively. 
In Eq. (l3T9l), 



(3.20) 



is the coupling constant of the phonon-mediated interac- 
tions between the electrons. Note that the combination 
Muj\_^ is the curvature of the interaction potential be- 
tween the carbon atoms, which has electrostatic origin; 
therefore, it does not depend on the carbon m ass M and 
scales as M(jj\^ 



Vag. Together with Eq. (2.18), this 



leads to the order-of-magnitude estimate 



h 



(3.21) 



Comparing with Eq. 
anisotropics (3.18) and 



, we see, that the bare 
due to short-range e-e and 
e-ph interactions are actually parametrically the same 
and can differ only numerically. 
The bare expressions 



3.16)-(3.19) determine the pa- 



rameters of the V = QHFM, provided one may neglect 
the effects of the interaction-induced electron transitions 
between different LLs, also known as "Landau level mix- 
ing" . As we find in the next section, this is not the case 
and, in fact, the parameters are drastically affected by 
the LL mixing effects. 
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IV. RENORMALIZATIONS OF PARAMETERS 
OF THE ly^O QHFM BY THE LANDAU LEVEL 
MIXING EFFECTS 

The weakness of the Coulomb interactions is generally 
a sufficient condition for the applicability of the QHFM 
theory. For graphene, this is formulated as e^/w ^ 1, in 
which case, since the typical Coulomb energy /Ib per 
orbital is much smaller than the LL spacing ei — \/2v/Ib, 
one could expect LL mixing effects to be inefficient and 
operate within the n — LL only. However, interac- 
tions are not weak in graphene: for suspended samples, 
e^/v sa 2.2 (taking v = lO^cm/s), which may be regarded 
as moderate strength. Besides, Coulomb interactions are 
known to be marginaP^ in graphene: they produce large 
logarithmic contributions in the diagrammatic series re- 
gardless of the their strength, even if weak. These loga- 
rithms come from the wide range of energies |e| ^ v/qq 
up to the bandwidth and involve many LLs. Therefore, 
taking the nonzero LLs into account is essential and the 
question arises whether the QHFM theory still holds for 
a realistic model of graphene, with e^/w ~ 1. 

In this respect, the large- iV ex pansio n in the number 
of "flavors" has gained popularit} !'^^ ! ^"^ ^ for graphene. 
This approach allows one to single out the leading in N 
diagrams, in each order in the bare Coulomb interactions, 
and perform partial summation of diagrams, formally 
analogous to the random-phase-approximation (RPA) se- 
ries. The physical justification of the method is that large 
N makes screening of the interactions especially efficient. 
This reduces the coupling constant from its bare value 
e^/v ^ 1 to that of the screened interactions ^ 1, 
resulting in the effectively weak-coupling theory. 

In this section, we will use the large-TV approach to 
systematically take into account the effects of LL mix- 
ing. In reality, iV = 4 in graphene due to two valleys and 
two projections of spin. Although this value is not par- 
ticularly large, one can still expect the large- approach 
to adequately describe the correlated physics in graphene 
for moderate Coulomb coupling e'^/v ^ 1. 

The effects produced by LL mixing can be identified al- 
ready in the second order in the interactions ( 2.9 1, ( 2.10 1, 
and (2.161. Fig.[5]shows the free-energy diagrams for the 
energy functional ( |3.11 ) , second order in the Coulomb in- 
teractions ( |2.9[ ), while Figs. |6] and [t] shows the diagrams, 
first order in the Coulomb and in either short-range e-e 
(2.10) or e-ph (2.16) interactions. As a starting point. 



one includes the contributions from all LLs in these di- 
agrams, using the full field operators (2.22). One then 
separates the contributions from n — and n ^ LLs 
in the electron Green's functions (solid lines), at which 
point the following observations can be made. 

The diagrams maintain the structure of the lowest- 
order diagrams in Fig. |4] with only n = LL present, 
but contain the blocks involving n 7^ LLs that repre- 
sent corrections to either (i) the Coulomb propagator - 
diagram in Fig. [sJaO); (ii) two-particle vertex functions 
- diagrams in Figs.[^(al), (a3)], |6[(bl), (b2), (b3)] and 





(al) 



(a3) 




FIG. 5: Diagrams for the energy functional (3.111, second 
order in the Coulomb interactions (wavy lines). Diagram (aO) 
represents the first-order large- A?^ correction to the Coulomb 
propagator. Diagrams (al)-(a4) diverge logarithmically, but 
cancel each other within the logarithmic accuracy. 



(bl) 



(b2) 



(b3) 



(b4) 







FIG. 6: Diagrams for the energy functional ( |3.11| l, first order 
in the Coulomb interactions (wavy line) and in either the 
short-range e-e or e-ph interactions (dashed line). Diagrams 
(bl)-(b4) diverge logarithmically and represent the lowest- 
order correction to the Hartree contribution [Fig.|4]^b)] to the 
anisotropy energy. 





(cl) 



(c3) 




{c2) 



{c4) 




FIG. 7: Same as in Fig. [6] but with respect to the Fock 
contribution [Fig. |4|c)] . 
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[7[(cl), (c2), (c3)]; (iii) or the self-energy of the n = 
LL Green's function - diagrams in Figs. 5ja4), |6|^b4), 
and[7Fc4). Without the magnetic field, the blocks (ii) 
and (lii) diverge logarithmically^, since Coulomb inter- 
actions are marginal; this remains true in the magnetic 
field. 

This pattern persists in higher orders, which allows one 
to formulate the general recipe for taking the LL mixing 
effects into account within the large- TV approach. One 
first performs a partial summation of the RPA-type series 
of blocks (i) . This produces the "dressed" propagator of 
the screened Coulomb interactions. One then performs 
the summation of the blocks (ii) and (iii). The lead- 
ing log-divergent diagrams can efficiently be summed up 
using the renormalization group (RG) procedure. One 
notices that, when combined together, the blocks (ii) 
and (iii) describe the critical ren ormal izations of the 
short-range e-e and e-ph interactionJ^^'^ by the screened 
Coulomb interactions. At the same time, the Coulomb 
interactions themselves are not renormalizecP^ e.g., the 
diagrams (al)-(a4) in Fig.jsjcancel each other within log- 
accuracy. 

With these steps performed, we find that, within the 
lar ge TV -approach, the v = QHFM theory [Eq. (|3.10[) 



standar d^^ l ^^^ l ^^ l expression 



'■B 
IGtt 



dqq^e — —V{uj^O,q) (4.2) 



in terms of an arbitrary potential. In Eq. (4.1 1, Il{uj,q) 
is the polarization operator of graphcne in the presence 



of the magnetic field. Since in Eq. (4.2) the frequency 
is constrained to a; = (the typical energy scales of the 
QHMF theory are <^ v/Ib) and the relevant momenta 
are g ~ 1//^ , one ha s to u se the exact expression for the 
polarization operatoiI^SESl. At filling factor v = zero 
temperature, and neglecting the minor corrections from 
the Zeeman effect, the polarization operator reads 



n(o,(z) 



N 



n'<0 



|ifn„'(q)|' 



(4.3) 



Here, i^„„'(q) [q = {qx,qy), q = |q|] are the graphene 
magnetic form-factors; they are expressed in terms of the 
conventional form-factors 



K„n' (q) 



+ CX3 



dxe"?-^ 



(3.151] does hold as a controlled approximation for 



graphene with moderate strength Coulomb interactions for the quadratic spectrum as 
e^/v ^ 1, but LL mixing effects result in crucial renor- 



(4.4) 



malizations of its parameters and u_l,z- These renor- 
malizations are considered in the remaining subsections. 

We emphasize that the separation of the contributions 
from n = and n 7^ LLs in the diagrams, with subse- 
quent classification of their blocks according to the types 
(i), (ii), (iii), is justified only in the weak-coupling large- 
N limit of the screened Coulomb interactions, when the 
real occupancy of the n ^ LLs is close to full or zero, 
1 or 0, for n < or 71 > 0, respec- 



i.e. 



tively. For example, for the general form of the second- 
order free-energy diagrams, with all LLs involved (left- 
most diagrams in Figs.[5j[6] and]?]), one cannot meaning- 
fully attribute each diagram to just one of the (i), (ii), 
(iii) classes. At the same time, virtual electron transi- 
tions between different LLs are quite efficient and result 
in the strong screening of the Coulomb interactions. 



A. Stiffness for screened Coulomb interactions 

The first consequence of the LL mixing effects is that 
the stiffness ps becomes suppressed due to screening, as 
compared to its bare value (3.161. The stiffness is ob- 

^ (a) , in which the 
Iq should be sub- 



Knn'iq) = -[i4:|„|,|„'|(q) -|-sgnnsgnn'^C|„|„ij„,|_i(q)], 

if both n ^ and n' ^ 0, and X„o(q) = ^|n|.o(q)- 

The polarization operator 11(0, q) depends only on the 
combination qls- For qls ^ 1, n(0,q) Nq/{16v) 
approaches its expression in the absence of the magnetic 
field. At arbitrary qlB, n(0,<?) can be calculated only 
numerically. 

Since the Coulomb potential Vo{q) has no scale and 
n(0, q) depends solely on qlB, the stiffness (4.2 1 scales as 



2nll 



Nl, 



-R{e^N/v) 



(4.5) 



where the dimensionless function Rje^ N /v) o f the cou- 
pling strength e^N/v is defined by Eqs. (4.2 1, (4.1 1, and 
(4.3). In the limit e^N/v <C 1 of negligible screening, 
y(0, q) ~ Vb((7), one obtains 



R{e^N/v < 1) 



TT e 

2~. 



tained from the Fock diagram in Fig. 
bare Coulomb potential Vo(<7) = 2ne 
stituted by the "dressed" propagator 



V{w,q) 



Vo{q) 



l + Vo{q)n{L,,q) 



(4.1) 



of the screened interaction (the second term of this series 
is shown in Fig. ^a.Oj). The stiffness is given by the 



recovering the expression (3.16) for the bare stiffness. 
The function R{e'^N/v) grows with increasing e^N/v 
and saturates to the maximum value R{oo) ^ 1 in the 
limit e^N/v ^ 1 of complete screening, when V{Q,q) w 

i/n(o,Q). 

We see that, upon taking the screening effects of LL 
mixing into account, the stiffness (4.5 ) retains its square- 
root scaling ps{B±)/ {2TTlg) cx \/B± with the magnetic 
field, but the numerical prefactor of the dependence be- 
comes suppressed. The main practical implication of this 
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concerns the activation transport through the bulk of 
the sample, since the gaps of the charge excitations are 
determined by the typical energy of symmetric interac- 
tions; e.g., the energy of the unit charge Skyrmions equals 



B. Renormalization of the anisotropy energies 

The second and a more physically significant ef- 
fect arising from nonzero LLs is the renormalization 
of the isospin anisotropy energies uj_^z [Eq- (3.14)] by 
the Coulomb interactions. 



The bare energies 



[Eqs. ( |3.17P , p.l8| , and p.l9[)] are determined by the 
diagrams (b) and (c) in Fig. |4j while the diagrams in 
Figs. |6] and [7] represent the lowest-order corrections to 
them. Since the latter diverge logarithmically, one is 
forced to sum up the whole series of log-divergent di- 
agrams. Within the logarithmic accuracy, this can be 
achieved by means of the RG procedure. It proves more 
efficient to consider the renormalizations of the bare cou- 
plings ga/B and f±z, rather than uf^^ directly, and then 
express the anisotropy energies m_l,z in terms of the renor- 

malized couplings. 

In the RG procedure for monolayer graphen^^^E^IlSl^ 



one starts with the Hamiltonian (2.1 1 of massless weakly 
interacting Dirac fermions. Integrating out the high- 
energy fermionic modes in the frequency-momentum 
space yields the RG eq uations for the involved cou- 
plings. One findd^^l^^Iini that the Coulomb interactions 
are marginally irrelevant and, consequently, weakly in- 
teracting electrons in monolayer graphene flow to a fixed 
point of noninteracting massless Dirac fermions. 

The RG analysis of the renormalizations of the short- 
range e-e and e-ph interactions by the Coulomb interac- 
tions was carried out in Refs. I44I45I Here, we recover 
the essential results and concentrate on the properties of 
key relevance to the = QHFM theory. Our notation 
differs from that of Refs. I44I45I 



1. Preparatory remarks 

Several comments are in order before we proceed. 

(a) In the perturbation theory diagrams, such as in 
Figs. [6]and[7j the large logarithms arise from the diver- 
gent sums over LLs e„ (if one first integrates over the 
frequency w), which have to be cut by the bandwidth at 
high energies, |eri| ;5 v/ao. Since such sums involve many 
LLs, the discreteness of the spectrum due to the magnetic 
field may be neglected: one may use the expressions for 
the Green's functions and polarization operators in the 
absence of the magnetic field, substituting the sums over 
LLs by the integrals over momenta q. At the lower limit, 
these integrals have to be cut by the inverse magnetic 
length q ^ ^/Ib, once the influence of the magnetic field 
becomes important. Hence, the arising logarithms are 



Therefore, the magnetic field does not affect the very 
structure of the RG equations of Refs. 1441451 yet defines 
a natural scale, at which the RG fiow stops. In the RG 
approach, the coupling constants 5a^(0 ^'^'^ f-Lz{l) ac- 
quire a dependence on the running length scale I. The 
RG fiow starts at the atomic scale I = a, where the cou- 
plings are equal to their bare values [Eqs. ( |2.10 1, (2.161, 
and (|3^], 



9aP 



(a) = gafs, f±zia) = f_ 



±z, 



and stops at the magnetic length I — Ib- The magnetic 
length Ib defines the scale, at which the renormalized 
anisotropy energies are to be determined. We define 
the atomic scale as a ~ oq, absorbing the ambiguity of 
the cutoffs in it. 

(b) Besides the renormalizations arising from the in- 
teractions in the process of "integrating out" higher en- 
ergy degrees of freedom ("mode elimination" part, in 
the terminology of Ref . [S3| , the full RG scheme also in- 
cludes the renormalizations due to rescaling of frequen- 
cies, momenta, and quantum fields (scaling, or "tree- 
level", renormalization). While the former represent ac- 
tual physical processes, described explicitly by the di- 
agrams (such as in Figs, [s] |6] and [?]), the latter are 
introduced in order to restore the original phase space 
and thus make comparison of the theories with different 
bandwidths meaningful. An important question arises, 
whether the latter, tree-level, renormalizations also have 
to be taken into account when determining the renor- 
malized anisotropy energies uj^ ^. Our understanding 
is the answer is negative, as they do not correspond to 
any physical processes. The following arguments can be 
given. 

(b.l) Suppose the logarithmic contributions due to 
the Coulomb interactions were absent or negligible: for 
the sake of argument, one may certainly consider a 
model with well-behaved finite-range symmetric interac- 
tions (e.g., Gaussian potential), or the Coulomb interac- 
tions so weak or the number of fiavors N so large that 
min(e^/?;, 1/A^) ln(ZB/ao) ^ 1. Then, all higher-order 
contributions to the anisotropy energy would be small 
compared to the first-order contribution, which is given 
by the diagrams (b) and (c) in Fig. |4] (the ultraviolet- 
divergent contributions, mentioned in Sec. |IIB[ from the 
short-range e-e interactions may also be assumed small, 
9apl{vao) ^ !)■ The renormalized anisotropy energies 
ua^^z would then be determined just by the bare expres- 
sions dSTTl), (Islsl), and (Isll), 



u^,z ^ u^^\. (4.6) 

At the same time, the full RG procedure would con- 
sist just of the tree-level (TL) renormalization and the 
short-range e-e and e-ph couplings would renormalize as 
.9a^(0 = 9a0j and fj^{l) = f±zj, since their scaling 
dimension is —1 (see paragraph (c) below regarding e- 
ph coupling). Therefore, taking the scaling renormal- 
ization into account, i.e., using the expression (3.181 for 
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the renormalized energies u± 

9al3 

bare energies u^^\ by the factor a/Is, 



in terms of the couphngs 



'^^{Ib) and /J^ (/s), one would have to multiply the 



a (0) 



This would, in apparent disagreement with Eq. (4.6 1, 



both drastically decrease the magnitude of the anisotropy 
energies and alter their dependence on the magnetic field. 

(b.2) Alternatively, such question and the extra fac- 
tor a/Is in u±,z never arise, if one starts with the "poor 
man's" formulation of the problem, as to sum up the com- 
plicated parquet series of log-divergent diagrams (Figs. [6] 
and[7]and all higher orders). One can then invoke the RG 
procedure, or rather, just its "mode elimination" part, as 
the latter is known to be an elegant way of accomplishing 
this task. 

We, therefore, conclude that the renormalized 
anisotropy energies u_l,z must be expressed through the 
couplings 5c(/3(^_b) and f±z{lB), in which only the renor- 
malizations arising from the "mode elimination" part of 
the RG procedure are taken into account: 



,(c-ph) 



(4.7) 



,(e-e) 



[gaoih) + 9az{h)]i ^=-L^Z (4-8) 



(o-ph) _ I±z{Ib) (c-ph) _ n 

L I — ,1-) • LL ^ — \j . 



(4.9) 



The couplings gapiP) and f±zil) differ from the couplings 
Qapil) and f±z{l) of the full RG scheme by the tree-level 
renormalization. 



5a/3(0 = 5a/}(0y> f±zil) = f±zil) 



(c) Finally, we comment on the peculiarities of e-ph 
interactions in graphene. Since the phonon dynamics is 
characterized by the phonon frequencies [fi — E2, Ai, 



Eq. (2.17l], the properties of e-ph interactions depend on 
the energy scale e at which considered. At e <C a;^, the 
phonon-mediated electron interactions are effectively in- 
stantaneous and, being also short-ranged, are quite anal- 
ogous to the short-range e-e interactions; in particular, 
they are irrelevant in the RG sense with the scaling di- 
mension — 1. At energies e S> o;^, retardation is strong 
and phonon-mediated interactions become marginaP^* 
they produce log-divergencies, which have to be cut by a;^ 
at the lower limit, yielding the logarithms ln[w/(aQa;^)]. 
As a result, in principle, e-ph interactions renormalize 
both themselves and the short-range e-e interactions. In 
practice, however, for the typical values of parameters in 
graphene, these renormalizations turn out to be numeri- 
cally smaller that those due to the Coulomb interactions, 
as can be inferred from the analysis of Ref. [45] There- 
fore, here, we neglect the renormalizations due to e-ph 



interactions. On the one hand, since e-ph interactions 
couple different valley channels, including these renor- 
malizations would significantly complicate the RG equa- 
tions. On the other hand, this should not alter the main 
conclusions of the RG analysis, which imply that essen- 
tially any algebraic possibility (i.e., including signs) for 
the isospin anisotropy energies {u±,Uz) could be realized 
in graphene. 



2. RG analysis 



We can now proceed with the RG analysi^^^'^. The 
one-loop renormalizations of the short-range e-e and e- 
ph couplings by the screened Coulomb interactions are 
actually described by the diagrams in Figs.|6]and[7l upon 
"dressing" the bare Coulomb lines. The corresponding 
RG equations for the coupling constants ^q,^ (0 reacP 



dffao(0 



— for a — x,y, z, 



(4.10) 



and 



dffaj_(0 



2F{wm9o.±il) - gcziOl (4.11) 



dgazjl) 
dC 



^ 4F{wmgo.z{l) ~ gc.±{l)], 



(4.12) 



for a — 0,x, y, z. Here, ^ = Inl/a. 
One observes that the couplings 

gao{l) = gao, a = x,y,z, (4.13) 

are not renormalized, whereas gax{l) = 5qi/(0 = .9qJ-(0 
and (jazil), Oi = 0,x,y,z, are. Different valley channels 
a are not coupled, but the _L and z sublattice channels 
do couple to each other within each valley channel. The 
reason for this are the properties of the electron Green's 
function: while a unit matrix in the KK' valley space, 
it has nontrivial matrix structure in the AB sublattice 
space. The symmetric Coulomb interactions themselves 
couple neither valley nor sublattice channels. 
In Eqs. (|4.11|) and (|4.12|). 



^, , " / , TT arccosw 
r (w) — — — 1 



(4.14) 



is a function of the dimensionless coupling constant 

ttN e2 



wil) 



8 v{l) 



(4.15) 



of the C oulomb interactions, plotted in Fig. |8] T he RG 
equatiorP^EsMl for is also determined by Eq. (4.141, 



dwjl) 



'F{w{l))w{l). 



(4.16) 
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F{w) ^ — 






w 



FIG. 8: The function F{w) [Eq. ( |47L4| )] of the Coulomb cou- 
pling constant w [Eq. ( |4.15 1] describes the renormalizations 
of the Dirac velocity v{l), w{l) itself, and the short-range e-e 
and e-ph couplings. 



w{l) 



l/a 



FIG. 9: Renor malization of the Coulomb coupling constant 
w{l) [Eq. (4.15l]. The bare value w(a) = 3.4 was used. 



The renormalization of comes from the renormal- 
ization of the Dirac velocity v{l), while the charge is 
not renormalized (the diagrams (al)-(a4) in Fig.[5]cancel 
each other within the log-accuracy). 



Equation (4.16) can be solved analytically only the 
limits of weak {w{l) ^ 1, interactions are not screened, 
F{w) ~ ) or strong {w{l) ^ 1, interactions are fully 
screened, F{w) ~ F{oo) = :^^) coupl ing. I n Fig. M we 



plot the numerical solution w{l) to Eq. (4.16 1 for thebare 
value w{a) = w w 3.4 obtained at v{a) — v ^ lO^cm/s, 
which should be typical for suspended graphene. As a 
general property, w{l) decreases [the velocity v{l) grows] 
with increasing the length scale /; therefore. Coulomb in- 
teractions flow towards weak coupling. At Ib/ci — 40 
(taking Ib ~ lOnm at B± « lOT and a — qq f» 2.5 A), 
we obtain w{Ib) ~ 2.1, i.e., for experimentally relevant 
values. Coulomb interactions remain in the intermediate, 
moderate strength, regime. 

Since the renormalization of the Coulomb coupling 
w{l) [Eq. (4.16)] does not involve the short-range e-e or 



e-ph interactions (we neglect the later, see paragraph (c) 
in Sec. |IVB 1| and Rcf. 45), one may regard F{w{l)) in 



Eqs. (4.11 1 and (4.12) as a known function of I. The sys- 



tem of the coupled equations (4.11) and (4.12) can then 
easily be solved by making a linear transformation to the 
coupling^^ 



-(0= o[5a.(0-5a±(0], 







\\\ 

\ \ 


y\ \\ \ 


\\'\ 









FIG. 10: RG flows [Eqs. ^4~17\ , ( 4~1^ , and ^4~19\ ] of the cou- 
pling constants §^±{1) and Qazi}) of the short-range electron- 
electron interactions (2.101. 



which translates the equations to 

= 0, — ^ 6F{w[l))ga-{l)- 



The solution of these equations is straightforward and, 
coming back to ga±{l) and gaz{l)^ one obtains 

ga±{l) = 7^{gaz +2ga±) - T^igaz - ga±)J^c-c{l / a, w) , 

(4.17) 

1 2 

gaz{l) = 7^{gaz + 2ga±) + T^igaz " ga±)^c-c{l / a, w) , 

(4.18) 

(4.19) 



Tc-dl/a, w) = exp 



' dl' 
-yGFiwin) 



for a — 0,x,y,z. The function J-c_(.(l/a,w) determines 
the strength of renormalization. It depends on the bare 
Coulomb coupling w as parameter and reaches its maxi- 
mum 



J--r(Va)=-^c-c(V«,oo)- 



in the strong coupling regime w ^ oo, when the interac- 
tions are fully screened. This provides an upper bound 
for the strength of renormalization. 

The RG flows in the space of couplings {ga±{l), gaz{l)), 
defined by Eqs. (14. 171) and (|4.18|), are plotted in Fig.flO 



Their key properties are as follows. Starting from the 
point {ga±ia),gaz{a)) = iga±,gaz) at the lattice scale 
I = a, the RG flow {ga±{l),gaz{l)) follows the straight 
line 



'^ga±il) + gaz{l) = 25„j_ + ga 



(4.20) 



This constraint implies that the sum of couplings 
I]/3=a;,y,z5a/3(0 = 2gaj. (0 + (0 different sublattice 
channels /3 is conserved under renormalization. 
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Further, depending on the relation between the bare 
values ga± and gaz, eventually, for strong enough renor- 
malization, the flow line ends up either in the quadrant 



{ga±il) > 0, gaz{l) < 0), if ga± > gaz 
or in the quadrant 



(ga_L(0 <0, >0), ifffa 



(4.21) 



(4.22) 



light red regions in Fig. [TOj That is, the initially al- 
gebraically larger (smaller) coupling always becomes or 
stays positive (negative). Deeper reasons for this inter- 
esting behavior could be sought in the chiral nature of 
the Dirac quasiparticles, which is the ultimate cause of 
coupling between _L and z sublattice channels. 

What concerns the renormalizations of e-ph interac- 
tions by the Coulomb interactions, the RG equations for 

{E2 phonons) and f±z{l) 



Ml, 



the couplings fz±il) 
{Ai , Bi phonons) reacP^ 



E2 



d/.±(/) 



- 2F{w{l))fzAl), 



AF{w{l))hz{l). 



(4.23) 



For e-ph interactions, different sublattice channels do not 
couple in the renormalization process. The solution to 
Eq. (|423l) is 



f±z{l) = f±zJ^o-phil/a,'i 



(4.24) 



case, when the bare values are close to the unstable line 
5aJ_ = gaz) to obtain 



,(c-c) 



1 2 



2ttII 3 



(gaz - ga±)J^c-c{lB/a,w), a^J-,Z. 



(4.26) 

From Eqs. ( |4.9[ ), (4.24), and ( 4.25[ ), the anisotropy ener- 
gies due to e-ph interactions equal 

^(o-ph) ^ (;^/^^ ^(c-ph) ^ Q (4 27) 

Final ly, w e also mention that the Zeeman energy ez 
[Eq. (3.13)] is not renormalizecP^ by the Coulomb in- 



teractions, since the Zeeman splitting term in Eq. (2.7) 
is scalar in KK' (g) AB space. 



C. Consequences of renormalizations of the 
anisotropy energies ?tx ^ for the u — QHFM 

The properties of the renormalized isospin anisotropy 



energies, obtained in Sec. IV B have very important con- 
sequences for the physics of the 1^ = QHFM. The con- 
clusions below constitute some of the key results of the 
present work. 



1. Magnitude of the anisotropy energies ux,^ 

The first consequence concerns the absolute values of 
the anisotropy energies u±^z- As follows from Eqs. (2.15 1, 
dSlSl), (ISlgl), and (E2TI) 



can be roughly estimated as 



the bare anisotropy energies 



J^c-ph{l/a,'w) = exp 



dl 



(4.25) 



Comparing Eqs. (4.191 and (4.25), we see that the 



renormalization of the e-ph coupling is parametrically 
weaker than that of the short-range e-e couplings: 
J^e-phil/o-Tw) = [J-c-cil/a,w)]'^^^ due to the factor 4 in- 
stead of 6 in the exponential. The maximum, as a func- 
tion of w, equals 



•^e'"p?(Va)=-^c-ph(Va,oo) 



Concluding this section, the renormalizations of the 
short-range e-e and e-ph couplings are described by 
Eqs. ( |413l ), ( |4T7t , dlts] ), ( |4l9| , (|X24]), and ( |425| . 
Since (jaoi^B) — gao, ct =_L,z, are not renormalized 
and, according to Eq. (4.8), gQp{lB)^ P =-L,2, do not 



affect the anisotropics, the main contribution to the 
anisotropy energies Ua °^ due to short-range e-e inter- 
actions comes from the couplings gaziln)- In the regime 
J'c-dlB/o-jw) ^ 1 of strong renormalization, one may 
neglect the first term in Eq. (4.18 1 (except for the special 



(o-c,0) 



(c-ph,0) 



Op 

ao \Ib , 

(4.28) 

They scale linearly with the perpendicular magnetic field 
B± and are on the same order as the Zeeman energy 
ez = ^J-bB ~ 0.7i3[T]K (for moderate tilt angles), as 
obtained earlier in Refs. 22 25 30 31, 
According to Eqs. ( |4.26[ ) and ( |4.27[ ), 



u 



(c-c) 



(Bi 



\u^rf\B^)\^,.,{lB/a,w) (4.29) 



and 



,(c-ph)/ 



(i?i) = uV"''''°^(5±)-^c-ph(;sA 



(4.30) 



i.e., the renormalized anisotropy energies are enhanced 
by the factors J-"o_o(Z_B/a, w) and J^c-ph(^s/a, it'), com- 
pared to the bare values and their functional dependence 
on B± is altered. 

we plot J^c- 



In Fig. 



11 



2(1 b/ a, w) and Fc-ph{lB / a,w) 

as functions of Ib jo, for the bare Coulomb strength w — 
3.4. Some values are [J^e.eils / a,w), J^c-pYi{Ib I oltw)) ~ 
(6.6,3.5), (19,7), (37,11) at Ig/a = 10,40,100, respec- 
tively. We see that, in a typical experimental situation 
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(c-c,0) 



Te-eilB/a, w) 



J^e-phils/a, w) 



Ib/cL: 



FIG. 11: The functions Ta-a{lB /cl,w) and J-e-ph(lB / cl,w) 



[Eqs. (4.191 and (4.25 1], characterizing the strength of renor- 
mahzations of the isospin anisotropy energies and u^^ '^"^^ 
[Eqs. (4.261 and (4.27l] due to short-range e-e and e-ph in- 



teractions. The value w = ?,A for the bare Coulomb cou- 



pling [Eq. (4.151] was used. In a typical experimental range 
Ib/u ~ 10 — 100, the anisotropy energies are enhanced by 



about one ord er compared to t heir b are values u, 



(g-g,0) 



and 



(o-ph,0) 



Eqs. (3.17l-(3.20l, and (4.28l] 



{Ib ~ lOnm and a « 2.5A gives Ib/cl « 40), the renor- 
malization of the anisotropy energy is expected to be very 
strong. Since the Zeeman energy ez = (J-bB is not renor- 
malized, u_l,z can easily exceed ez by one order. This 
implies that the isospin anisotropics play a major, per- 
haps, more significant role than the Zeeman effect in the 
physics of the ly ^ QHFM. 



2. Signs of the anisotropy energies ux,z 

The second consequence concerns the possible signs of 
the anisotropy energies u±^z- According to Eq. (4.26), in 
the regime of strong renormalization J-c-c{Ib/o., w) 3> 1, 
the sign of the anisotropy energy Ua~'^\ a =-L,z, is de- 
termined by the relation gaz ^ 5a_L between the bare 
couplings. This follows from the discussed p eculia r be- 
havior of the RG flows, see Fig. 10 and Eqs. (4.21) and 



(4.22). This property suggests that. 



essentially, any sign combination 



(c-c) „ 

u± ' ^ 0, 



i^r"^ ^ 



(4.31) 



5^/(20 



JaO 



> 0, a =_L, z, can also only be 
,(1) 



< gLs for e^/v < 1). 

(e-e,0) 



U, 

positive {gao 

If the bare expressions m^'^'"'' provided correct values 
for the anisotropy energies, this would significantly re- 
strict the variety of possible orders of the i/ = QHFM 
that could be realized in the system. Anticipating the 
results of Sec. |Vj the fully spin-polarized ferromagnetic 
phase would be the only phase favored by the short- 
range e-e interactions. However, the positiveness of the 
bare couplings g^^l > does not impose any sign restric- 



tions on the renormalized energies u 



(c-c). 
±,z ■ 



according to 

Eq. ( |4.26[ ), the signs of Ua'*''' {a =-L, z) are determined 

by the relations gaz ^ g^^l^ between the couplings, and 
not their signs, and there seems to be no physical restric- 
tion on the relative values ga] / g^^}^ . 

This demonstrates that, even if the signs of the bare 
couplings gap {a,/3 =_L,z) were fixed to positive by the 
repulsive nature of the Coulomb interactions, the renor- 
malized anisotropy energies can still come in any 
possible sign combination, which brings us to the above 
statement. We also do not expect the sign restriction on 
ga/s (q^, /3 =-L, z) to necessarily remain in the regime of 
stronger interactions, e'^/v ^ 1: the sign change of some 
couplings could occur already due to renormalizations at 
the lattice scale, see Sec. |II B[ 

What concerns the anisotropics due to electron inter- 
actions with in-plane Ai,Bi phonons, as follows from 
Eqs. (|3.19|) and (|4.27|), the energy 



P*^' < 



(4.32) 



always remains negative. The negative sign is a conse- 
quence of the attractive nature of the phonon-mediated 
interactions and appears to stay preserved in the renor- 
malization process. We also mention that electron inter- 
actions with the out-of-plane phonons result in a negative 
anisotropy energy ui'^'^^^ < 0, which is, however, small, 

,(o-ph)| 



,(o-ph)| 



u'z ' 'I ^ *"*'| and neglected here. 

The discussion of the implications of the properties 
( |4.31[ ) and ( [4.32^ for the physics of the = QHFM will 
be continued in Sec. IV Bl 



of the anisotropy energies due to short-range e-e interac- 
tions could be realized in a real graphene system, regard- 
less of the potential restrictions on the signs of the bare 
couplings gap, a, /? =-L, z. 

To elaborate on this statement, it is instructive to con- 
sider the weak-coupling limit e^/w <C 1. In this regime, 

the bare couplings gap — g^^p, a,/3 —J-,z, obtained in 
the first order in the Coulomb interaction, are given by 
Eq. (2.12). The expression (2.12) in the form of the elec- 
trostatic energy leads us to conclude [Eq. ( |2.13 1] that 



the couplings g^^^ > 0, a, l3 —J-,z, are positive, i.e., 
e-e interactions in these channels are repulsive. Then, 



according to Eq. (3.181, the bare anisotropy energies 



V. PHASE DIAGRAM FOR THE i/ = 
QUANTUM HALL FERROMAGNET 

In this section, we obtain the phase diagram of the 
h' — QHFM in graphene in the presence of the isospin 
anisotropy and Zeeman effect, as described by the en- 
ergy functional E[P{t,r)] [Eqs. ( |3.11[ )-( |3.15p ]. We will 
consider only the "classical" ground states, i.e., time- 
independent configurations P{r) that minimize the en- 
ergy functional; effects of thermal fluctuations are be- 
yond the scope of this paper. In the bulk of the sample, 
the functional -^[^(r)] is minimized by a spatially ho- 
mogeneous configuration P{r) = P that minimizes the 
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u = QHFM state 


single-particle spinors Xa.b 


order parameter P 


anisotropy energy £o{P) 


spin-polarized isospin-singlet 


\n) ® |s), 1 - n) ® |s) 




—2u± — Uz 


isospin-polarized spin-singlct 


|n) ® |s), |n) (g) 1 - s) 


Pn® i 


«_L(n^ -I- Uy) + UzU^ 


P"" state 


|n) (g) |s), 1 - n) ® 1 - s) 


Pn ® Ps + P-n ® P-s 


/ 2 1 2 \ 2 



TABLE I: Relevant states of the v = Q QHFM. For given values of {u±, Uz), the anisotropy energy fo(P) is minimized by one 
of these states, see Tab. [II] 



sum 



£{P) = £^{P)+£z{P) 



(5.1) 



of the anisotropy and Zecman energies. Which P de- 
livers the minimum of £{P) depends on the signs of the 
anisotropy energies uj__z and their relations between each 
other and the Zeeman energy ez ■ The results of the pre- 
vious section suggest that essentially any scenario for 
(u^,Uz) could take place in real graphene. Therefore, 
here, we consider all possibilities and obtain a generic 
phase diagram in the space of parameters (u_l, u^^tz)- 



Relevant states 



According to Sec. Ill A the order parameter matrices 
P of the V — Q QHFM in graphene form a represen- 
tation of the eight-dimensional Grassmanian manifold 
Gr(2,4). For the purpose of finding the ground states, 
however, considering the most general parametrization of 
P is not ne cessary, since the anisotropy terms [Eqs. (3.14) 
and (3.15)] act explicitly on the isospin only. Here, we 



will consider two simpler six-dimensional submanifolds of 
Gr(2,4) of physical relevance to the problem. We will find 
that all possible ground states belong to these subsets. 

First, consider a family of states [see Eqs. (3.5) and 
(3.7), direct products are in KK' ® s space], 



n„ 



Xb 



P" = Pn, ®Ps+ Pn, ® P- 



(5.2) 
(5.3) 



in which two electrons occupy the states Xa,b with oppo- 
site spin (±s) and arbitrary isospin (n^^b) polarizations. 
Here and below, 

Ps = |s>(s| = i(i + rs) 

is the density matrix of the spin (isospin, if in KK' space) 
in the state |s); ria^;,, s, and n, Sa^b below, are unit vectors 
defining the spin or isospin states. Calculation of the 
anisotropy energy [Eqs. (3.14) and (3.15)] of the state 
(5.3) gives 



£o{P'') = U^{naxnbx + nayUby) + U^Ua^Ubz (5.4) 

As an important special case of P" states, for coincid- 
ing isospins iIq = n?, = n, one obtains a fully isospin- 
polarized (IP) spin-singlet state 

piP = P„ «) i (5.5) 




FIG. 12: Charge-density-wave (CDW) phase of the u = 
QHFM. 




FIG. 13: Kekule-distortion (KD) phase of the u = QHFM. 



with the anisotropy energy 

£o{P^^) = uj_n]_ + Uzul, = nl + n^. (5.6) 

The observable order of such state depends on the orien- 
tation of the isospin n. In particular, when the isospin 
is at the poles of the Bloch sphere, n = in^, = 
(0,0, 1), the state has a charge-density- wave (CDW) or- 
der (Fig. 12), with both electrons per orbital occupying 
the same sublattice, and energy £^^^ = Uz- When the 
isospin is on the equator of the Bloch sphere, n = nj^ = 
(cos sin (^„, 0), the state has the Kekule-distortion 
(KD) order, represented schematically in Fig. 13 and 
energy £^^ — u±. The angle tfn is related to the orien- 
tation of the atom displacements in Fig. |3| right). 

The second important family of states is "dual" to P", 



Xa = |n) (g) |sq), Xfc = I - n) (g) S{,), (5.7) 



P' = P^®P,^+P_^®P,^^ (5.8) 

Here, two electrons occupy the states Xa,b with oppo- 
site isospin (±n) and arbitrary spin {sa^b) polarizations. 
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FIG. 14: Spin-polarized ferromagnetic (F) phase of the u = 
QHFM. 



FIG. 15: Antiferromagnetic (AF) phase of the u = QHFM. 



Calculation of the anisotropy energy gives 



{uj_n^+u,n^) . 

(5.9) 

As a special case of P'* states, for coinciding spins Sa ~ 
Sfc = s, one obtains a fully spin-polarized isospin-singlet 
state 



1® Ps 



(5.10) 



with the ferromagnetic (F) order (Fig. 14) and anisotropy 
energy 



-2u± — Uz 



(5.11) 



Note that in the isospin-singlet F state, each isospin chan- 
nel a = x^y,z contributes —Ua to the anisotropy en- 
ergy (13. 141). 



An intersection of the subsets (5.3) and (5.8) is the 
family of states 



Xa 



n 



pn 



|s), Xb = 



Pn «) Ps + P-r 



P-s 



(5.12) 



Here, two electrons have simultaneously opposite spin 
(±s) and isospin (±n) polarizations. For brevity, we will 
denote the states (5.12) with n = n_L and n = ±n^ 



as P"* and P"*, respectively. At n = in^, electrons 
with opposite spin polarizations ±s reside on different 



sublattices and P"** state has an antiferromagnetic (AF) 
order (Fig. 15). 



The state (5.12) has the anisotropy energy 



-u±n']_ 



(5.13) 



5.12) 



as can be obtained from either Eq. (5.4) or (5.9). Com- 
paring Eqs. (5.9), (5.11) and (5.13), we also notice 
that, for arbitrary orientations of spins and Sf,, the 
anisotropy energy (5.9) of the state (5.8) is a linear com- 
bination of the energies of the F and P"" states, with 
the weights determined by the product s^st,, i.e., by the 
angle between the spins. 

The properties of the states (5.5), (5.10), and 
are summarized in Table U 



B. Phase diagram neglecting the Zeeman effect 

To get a clear understanding of the orders favored by 
the isospin anisotropy, let us first neglect the Zeeman en- 
ergy £z{P) completely and find the states P that mini- 
mize the anisotropy energy £o{P) [Eqs. (3.14) and (3.15 )] 
alone. 

We accomplish this by noticing the following property. 
For a given isospin channel a = x,y, z, the function ta{P) 



Eq. (3.15)] belongs to the range 



-2 < ta{P) < 2. 

The maximum ta{P) = 2 is reached at the IP 
state (5.5) with n parallel to the a axis, ria = ±1- For 
the remaining components a, one has ta{P) — 0. 

The minimum ta{P) = — 2 is reached within the subset 



(5.8), when 



(l-n2)(s+s_-l)=0, 

i.e., in two cases. 

(i) In the first case, S-|_ = s_, and the minimum is 
reached at the F state (5.10) (hence, n can be arbitrary). 
For the remaining components a, one also has ta{P) = 
—2, see also the comment after Eq. ( |5.11 ). 

(ii) In the second = ±1. For the remaining 
components a, one has ta{P) — —(1 + s+s_), for which 
the maximum ta{P) = is reached at s_|_ = — s_, i.e, at 
the P„^ (a =_L,z) state, and the minimum ia(P) = —2 
- at the F state. The latter just brings us back to the 
case (i). 

Using these properties, we arrive at the conclusion 
that, for given signs and ratio of u±^ and Uz^ the 
anisotropy energy f <> (P) is minimized by one of the above 
states - F, IP, or P"'' - that either minimize or maximize 

ta(P). 

The four possible cases (u_l ^ 0, ^ 0) of sign combi- 
nations are considered below; three of them split into sub- 
cases, depending on the relative absolute value |m_l/uz|. 

(-f-H): > 0, Uz > 0. The anisotropy energy 

is minimized by the F state (5.10), which minimizes 
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sign combination 


u± > 0, Uz > 


itx > 0, < 


u± < 0, Uz > 


u± < Q, Uz < Q 


subcase 






Uj_ < \Uz\ 


\U±\ > Uz 


M± < Uz 


\u±\ > \Uz\ 


\ui_\ < \Uz\ 


state minimizing £o{P) 


F 


F 


CDW 


KD 


AF 


KD 


CDW 


minimal So{P) 


—2u± — Uz 


— 2m_l + \uz\ 


— \Uz\ 


- -"a 


— Uz 


-\ui_\ 


-\Uz\ 



TABLE II: Minimization of the anisotropy energy £<,{P) [Eqs. ( |3.14| and ( |3.15[ )] 




FIG. 16: Phase diagram of the v — QHFM states mini- 
mizing the isospin anisotropy energy £o{P) [Eqs. ( 3.14[ ) and 
(3.15 1], in the space of the anisotropy energies {u±,Uz). Phys- 
ical orders of the phases are shown in Figs. |12p5| 



ta{P) = ^2 simultaneously for all a = x,y,z, and gives 
= -2uj_ - u^. 

( ): u± < 0, Uz < 0. The anisotropy energy is 

minimized by one of the IP states (5.5), which maximize 
taiP) = 2. For > \uz\, £<,{P) is minimized by the 
KD state with — — whereas for < \uz\ - 
by the CDW state with ffDW ^ -\uz\. 

( — h): u_L < 0, Uz > 0. The anisotropy energy is 



KD 



minimized by either the KD state with 
which minimizes the £±{P) — hu±[tx{P) + ty{P)] part 
of or the AF state [Eq. ( |5.12p with n = ±n^] with 

= —Uz, which minimizes the £z{P) = \uztz{P) part 
of £o{P), £o{P) = £±{P) + £z{P)- Comparing these en- 
ergies, we obtain that for > Uz, the KD state is 
realized, whereas for |u_l| < u^, the AF state is realized. 
The F state, which also minimizes £z{P), cannot be re- 
alized, since £^ = 2\u±_ \ — Uz > £^^ in this regime. 

(H — ): u_L > 0, < 0. In this case, one has to com- 
pare the energies of three states: the CDW state with 
^■CDW _ — which minimizes the £z{P); the F state 
with £^ — —2u± + \uz\, which minimizes £±{P) but 
maximizes £z{P) at the same time; and the P"* state 
[Eq. (|5.12|) with riz = 0] with £<,{Pl') = which 



has a higher £±{P) than the F state, but zero £z{P), on 
the other hand. Comparing these three energies, we find 
that only the first two states are realized: for u± > \uz\, 
£o{P) is minimized by the F state, whereas for u± < \uz\ 
- by the CDW state. The P"* state is not realized since 
5o(Pr') > £o for u_L > \uz\ and 5o(Pr) < ^o^"^ ^r 
u± < \Uz\. 



state 


order parameter P 


S4P) 


symmetry 


F 




—2u± — Uz 


SU(2), 


KD 


Pn_L ® i 


u± 




CDW 


P±n, ® i 


Uz 


Z2 KK' 


AF 


Pn, ® Ps + P-n, ® P-s 


-Uz 


SU(2), 



TABLE III: States minimizing the isospin anisotropy energy 
£o{P) [Eqs. (3.141 and (3.15l] and forming the phase diagram 
in Fig. 16 The last column (symmetry) denotes the symme- 
tries of the ground states in the isospin {KK') and spin (s) 
spaces. 



These cases are summarized in Tab. [TTj Together, they 
combine into the phase diagram of the states of the v = 
QHFM that minimize the isospin anisotropy energy, 
plotted in Fig. |16[ The phase diagram consists of four 
states with the following orders: F, AF, CDW, and KD, 
T2][15 



shown in Figs. 



The phases are separated by the 



first-order transition lines: 



> - between the F and CDW phases; 



u± — Uz < - between the CDW and KD phases; 



< - between the KD and AF phases; 



u± = 0, Uz > - between the AF and F phases. 

These transition lines come together at the origin 
{u±,Uz) = (0,0), where the anisotropy energy vanishes 
and all orders P have the same energy. 

Without the Zeeman effect, the phases have the fol- 
lowing symmetries: the F and AF phase are SU(2)- 
symmetric in the spin space [the spin polarization s in 
Eq. (5.10) or (5.12| can be arbitrary]; the KD phase has 
U(l) symmetry [n — n± = (cos 93„, sin 93„, 0) can have 
arbitrary orientation ipn in the Uz = plane], and the 
CDW phase has Z2 symmetry (n ~ ztUz, according to 
the occupation of either A ot B sublattice). Exactly at 
the phase transition lines, the symmetry of the ground 
state becomes higher; we do not attempt to study the 
details of the transitions here. The key properties of the 
phases are summarized in Tab. |1II[ 

Let us come back to the microscopic origins of 
the i sospin anisotropy. The conclusions of Sec. IV C 
[Eq. (4.31)] imply that any possibility for the signs and 



relative values of the anisotropy energies u^^ could be 
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realized and, therefore, any phase on the diagram in 



Fig. 16 could be favored by the short-range e-e interac- 
tions alone. The reason for that are peculiar properties 
of the renormalizations of the short-range e-e interac- 
tions in graphene, which allow for sign changes of the 
coupling constants (Fig. 10 1, switching the interactions 
from repulsive > 0) to attractive < 0) in cer- 
tain channels. 



In contrast, according to Eq. (4.32), the leading e-ph 
interactions always favor the Keku le dis tortion order, in 
agreement with earlier prediction^^oEI] Unlike short- 
range e-e interactions, different sublattice channels do 
not couple in the renormalization process and e-ph cou- 
plings retain their negative sign, characteristic of attrac- 
tive interactions. We also mention that e-ph interactions 
with the out-of-plane phonons, neglected here as weak, 
result in ui'^'^^^ < and favor CDW order, in line with 
Ref.El 

Of course, the values of the bare couplings ga/3 are 
determined by the details of e-e interactions and band 
structure at atomic scale and should be a robust mate- 
rial property. Therefore, in real graphene, one can ex- 
pect one particular situation for the anisotropy energies 



("1 



(e-e) „,(e-c) 



) to be realized [which, in the experimen- 



tally relevant regime of strong renormalizations, cannot 
be changed by varying B±, see Eq. ( |4.26 1 and discus- 
sion in Sec. VE and one certain order to be favored. In 



this sense, a more accurate formulation of the statement 



(4.31) is that one cannot theoretically rule out any pos- 
sibility for u^'^'', based just on the repulsive nature of 
the underlying Coulomb interactions. Doing so requires 
a reliable numerical estimate for the bare couplings. 

In the absence of such an estimate, we point out one 
case that may be the most relevant to the real sys- 
tem. For t hat, we turn again to the first-order expres- 
sions (2.12) for the couplings g^^j, a,/3 =_L,z. There, 
Pap(y), oi-tP — x,y, and Pzzir) are the staggered-type 
densities given by the linear combinations of the prod- 
ucts u^Airju^'Air) and u^B{r)u^,'B{r), P' = K, K' , of the 
Bloch wave-functions. On the other hand, the densities 
Paz(f) and pzair), OL — x,y, are determined by the over- 
laps Ufj_A(f)Ufj_iKB{r) of the Bloch wave- functions peaked 
at the atomic sites of different sublattices. Therefore, it 
would be reasonable to expect the couplings g'^^l and 
to be smaller than gi]) and g^^]_, i.e.. 



iJzz I i/_L_L ^ y±z' yz± ^ 

(This certainly becomes true in the limit of "strongly 
localized" atomic orbitals, when the size of the atomic 
wave-function is much smaller than the bond length.) In 
this case, according to Sec. IV B gzzils) > will stay 
positive and grow upon renormalization, while g±z{lB) 
will turn negative and grow in absolute value. This re- 
sults in the signs 



Ji) „(i) 



uf'^ < 0, 



,(o-o) 



Ml 

Ml 



J 




Sb x:^- — i 
/ 








\x 

\ ^^^^^ 









FIG. 17: Canted antiferromagnetic (CAP) phase of the = 
QHFM. 



of the anisotropy energies [Eq. (4.8)]. The same holds for 
the total anisotropy energies, 

u± < 0, Uz > 0, 



once e-ph interactions [Eq. (4.32)] are included. In this 
case, according to Fig.[T6]only either KD or AF phase can 
be favored by the isospin anisotropy (the latter becomes 
canted in the presence of the Zeenian effect, see Sec. V C I. 



C. Phase diagram in the presence of the Zeeman 
effect 

We now take the Zeeman effect into account and find 
the ground states P that minimize the the sum £{P) 
[Eq. (5.1 1] of the anisotropy and Zeeman energies. The 



question is how the phase diagram in Fig. [I6| is modified 
by the Zeeman effect. 

We first note that the spin-singlet CDW and KD 
phases are unaffected by the Zeeman field and their total 
energy £{P) is equal to the anisotropy energy, 



CDW 



KD 



(5.14) 



One the other hand, the "spin-active" F and AF phases 
are affected by the Zeeman field and, therefore, their 
whole sector has to be reconsidered. Since both of these 



states belong to the family (5.8), it is sufficient to mini- 
mize the energy 

£{P'') = -Uj_ -U^- U±SaSb - ez{Saz + S^z) (5.15) 

of P" state with n = but generally noncoUinear spins 
Sa^b- One can easily check that for a given angle 26s 
between the spins, SaS^ — cos(26's), Q < 9s < 7r/2, which 
fixes the anisotropy energy £o(^*), the Zeeman energy 
£z{P^) is minimized by the spin orientations 



Sa,6 



si cos(^s,±^l - si sin (/5s, 



(5.16) 



Sz — COS 6s, that have equal projections on the direction 
of the magnetic field (z) and are antiparallcl in the per- 
pendicular (xy) plane, as shown in Fig. 17 The total 
energy (5.15) then equals 



> 



£isz) 



u^{2sl 



1) 



2ezSz 



(5.17) 
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Minimizing it with with respect to Sz, we obtain that, 
for u± < —ez/2, a canted antiferromagnetic (CAF) state 
with the optimal angle 6** between the spins (Fig. 18), 



and energy 



7CAF 



2\u^\ 



(5.18) 



(5.19) 



is realized, whereas for uj_ > —6^/2, the fully spin- 
polarized F state (0* = 0, s* = 1) with the energy 



-2m I 



2ez 



(5.20) 



is realized. 

We see that in the presence of the Zeeman effect, the 
total energy (5.191 of the CAF state is indeed smaller 
of the AF state with antiparallel 
spms: by tormmg a noncoUinear orientation (Fig. 17), 



( |5.19 ) 

than that £'=^^^^^ 

by forming a noncoUinear orientation (Fig 
electrons lose some of the anisotropy energy, but gain 
more in the Zeeman energy. Therefore, the AF phase in 
Fig. 16 is completely substituted by the CAF phase, in 

depends on the ratio €z/\u±\. 
7t/2, si = 0) is reached 



which the optimal angle 0* 
The antiparallel orientation [9] 
asymptotically for ^ ez- 

We can now determine the boundaries between differ- 
ent phases. As obtained above, the F and CAF phases 
are separated by the line 



2 ' 



(5.21) 



The separation line between the CDW and KD phases 



Eq. (5.14)] remains at 



Next, comparing the energy (5.19) of the CAF state with 



that (5.14) of the KD state, we obtain the separation line 



-ez/2 




FIG. 18: Phase diagram of the v = Q quantum Hall ferro- 
magnet in monolayer graphene in the presence of the isospin 
anisotropy and Zeeman effect, in the space of the anisotropy 
energies {u±,Uz) and with the Zeeman energy ez as parame- 
ter. Ph ysic al orders of the phases are shown in Figs. |12| |13| 



14 



and 



17 Top graph shows the optimal value s* [Eq. ( 5.18 1 
of the spin projection on the direction of the magnetic field 
in the CAF and F phases. 



state 


order parameter P 




symm. 


F 


i (8) Ps.- 


— 2u±— Uz — 2€z 


none 


KD 




u± 




CDW 






Z2 KK' 


CAF 




^ 21^^x1 


U(l). 



TABLE IV: States minimizing the sum £o(P) + £z{P) of 
the isospin anisotropy and Zeeman energies and forming the 
The last column (symm. 



(5.22) phase diagram in Fig. 18 



denotes 

the symmetries of the ground states in the isospin (KK') and 
spin (s) spaces. 



(5.23) 



Analogously, the phase boundary between the F and 
CDW phases is now given by 



Uj_ 



Uz = -ez- 



(5.24) 



All four lines (5.21 )-( 5.24) come together and terminate 
at point (uj_, Uz) = — (1, 1)6^/2. 

This forms the phase diagram of the v = quantum 
Hall ferromagnet in monolayer graphene in the presence 
of the generic isospin anisotropy and Zeeman effect, plot- 
ted in Fig. [18] This diagram constitutes the key result of 
the present work. The total energy (5.1 ) is minimized by 



the states with one of the following orders: spin ferromag- 
netic (F), charge-density- wave (CDW), Kekule distortion 
(KD), or canted antiferromagnetic (CAF). As compared 
to the situation without it (Fig. 16), the Zeeman effect 



(i) substitutes the AF phase with antiparallel spins by 
the CAF phase with noncoUinear spins and (ii) naturally 
widens up the region of the F phase in the {u±,Uz) plane. 
The discussion of the microscopic origins of the isospin 
anisotropies done in Sec. |VB] can be directly carried over 
here. The key properties of the phases are summarized 
in Tab. JVl 

As far as the symmetries are concerned, the fully 
isospin-polarized CDW and KD phases, unaffected by 
the Zeeman field, retain their U(l) and Z2 degeneracies 
of the isospin orientation. At the same time, the F phase 
becomes nondegenerate, with the spin s = Sz = (0,0, 1) 
directed along the Zeeman field (i.e., the total magnetic 
field), and the CAF phase is U(l)-degenerate with re- 
spect to simultaneous rotations of the spins Sa^b about 



the direction of the Zeeman field [angle (ps in Eq. (5.16)]. 
The continuous U(l) degeneracies of the CAF and KD 
phases could be subject to thermal fluctuations, which 
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we do not addressed here. 

Note that, while all phase transitions in Fig. [16] and 
the rest of the transitions in Fig. 18 are first-order (black 
lines), the CAF-F transition [Eq. (5.21 ), dashed blue line 
is second-order: 



in Fig. 18 



upon increasing ezl\u±\ 



the 



CAF phase continuously crosses over to the F phase, 
as the AF component (~ ^1 — s*^) of the CAF order 
parameter gradually decreases, while its F component 
s*) grows; eventually, at the CAF-F transition line 
(s* = 1), the AF component turns to zero, while the F 
component saturates and experiences a jump in deriva- 
tive (see top graph in Fig. 18). 



D. Relation to earlier works 

Here, we discuss the connection of our results to earlier 
related studies. 

In Refs. 1221291 the lattice effects of e-e interactions on 
the V = Q QHFM were studied using the tight-binding 
extended Hubbard model with adjustable interactions 
at the lattice scale and asymptotically Coulomb inter- 
actions at large scales. Within this model, the overlap 
of the orbitals at different atomic sites is exactly zero 
and, neglecting renormalizations, the anisotropy energy 



,(<=-c) 



vanishes exactly and only u'z "'^ is present. The 

< are realized when, roughly, 
the sum of three nearest-neighbor repulsions is smaller 
or greater than the on-site repulsion, respectively. Ac- 
cordingly, in Ref. [211 in the presence of the Zeeman ef- 
fect, the competition between the F and CDW phases 
was predicted, which agrees with the phase diagram in 
Fig. [18] at the line uj_ = 0; the transition point is given 
. In Ref. Uni the comparison between the 
AF phases was done using numerical mean- 



cases Uz > and Uz 



by Eq. ( |5.24 

CDW, F, and 



field analysis. The conclusion was reached that, ignoring 
the Zeeman effect, depending on the details of interac- 
tions at lattice scale, either CDW or AF phase is fa- 
vored, while the F phase has higher energy. This is also 
consistent with the phase diagram in Fig. [16] at u±^ = 0. 
Note that, at u_l = and Uz > 0, the system is right 
at the transition line between the AF and F phases and 
therefore, even minor perturbations (numerical or other) 
would drive the system into one of the phases. 

One may also draw certain parallels between the v = 
QHFM in graphcne and the semiconductor quantum Hall 
bilayers"*^ ^^^'^ (QHB) at the total filhng factor i'qhb = 
2. There, the role of the isospin is played by the layer de- 
grees of freedom. The leading anisotropy comes from the 
difference between the Coulomb interactions within and 
between the layers. The resulting "capacitance" effect is 
described by the anisotropy Uz > 0, while a minor u± > 
due to the Coulomb interactions is usually neglected in 
theoretical studies. 

The proximity of the layers results a finite overlap of 
the different-layer wave- functions and a possibility of tun- 
neling. In the QHFM theory, this is described by an ex- 
tra isospin "Zeeman" term £t{P) = —etti[TxP], e* > 0, in 



the energy £{P) [Eq. ( jS.l] )], i.e., tunneling by itself favors 
the isospin-polarized (IPx) state P^p=" = P^^ (g) 1* with 
the isospin n^: — (1,0,0) along the x direction. Remark- 
ably, in addition to the F and IP^; phases favored by the 
Zeeman effect and interlayer tunneling, the CAF phase 
with the spin polarizations of the layers as in Fig. [17] was 
also predictecpS'^ to exist in a finite-size region of the 
phase diagram between the F and IPx phases. According 
to Refs. [55H571 the antiferromagnetic coupling between 
the spins in different layers favoring the CAF phase has 
a super-exchange nature and arises from the correlated 
two-particle tunneling processes. 

It should be emphasized that the physical origin of KD 
and CAF phases in graphene and IP^ and CAF phases in 
QHB is different. On the one hand, the possibility of tun- 
neling between the layers is a necessary condition for the 
existence of both IP^; and CAF phases in QHB. An anal- 
ogous ingredient is absent in the bulk of real graphene 
samples: creating the isospin "Zeeman" field directed in 
the xy isospin plane requires a static "nano-engineered" 
Kekule distortion. Fig. jsjright). On the other hand, in 
graphene, the KD phase is favored by large enough neg- 
ative anisotropy uj_ < 0, > while the AF phase 
(CAF, in the presence of the Zeeman effect) is favored 
by the combined anisotropics Uz > and u± < at 
Uz > —u±. The negative anisotropy uj_ < arises from 
the attractive interactions in _L-isospin channel, provided 
by either the e-ph interactions or short-range e-e interac- 
tions that turned attractive upon renormalization. Such 
mechanisms of u± < are minor or absent in QHB. 



E. Inducing phase transitions in the i/ = QHFM 

Turning to potential practical applications of the 
present theory, two important questions can be ad- 



dressed: (i) which phase in Fig. 18 is realized in the ex- 
perimentally observed insulating 1^ = state; (ii) whether 
one can induce transitions between different phases in a 
real system. We address the latter question in this sub- 
section and the former in Sec. I VII 

The ground state of the v = QHFM is determined 
by the relations between the parameters {u±,Uz,ez), the 
anisotropy u±^z and Zeeman ez energies, and the ques- 
tion is whether one can change these relations in the ex- 
periment. Changing the bare couplings g^p a, l3 =_L, z, 
in order to modify the relation between u± and m^, seems 
quite challenging as they are determined by the details 
of the interactions and band structure at atomic scale. 
Varying the magnetic field, in magnitude or orientation, 
is then virtually the only practical option. 

Let us first neglect the Zeeman effect and inquire if 
the relation between u±{B±) and Uz{B±) could be mod- 
ified by varying the magnitude of the perpendicular mag- 



netic field B±. As the RG analysis of Sec. IV B shows, 
the s hort- range e-e couplings gagils) [Eqs. (4.171, (4.181, 
and (4.19 )] do change their signs and relative values upon 



renormalization, and, as a result, so do the anisotropy 
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FIG. 19: Phase transitions induced by tilting the magnetic 
field, which increases the Zeeman energy ez relative to the 
anisotropy energies ux,z. The transitions from the CAF 
(point a) or CDW (point c) to the F phase occur directly, 
while the transition from the KD (point b) to the F phase can 
occur only through the CAF phase. 



energies u±^z{B±) [Eqs. (4.7) and (4.8 1] upon varying 
the magnetic field B±. However, according to the es- 
timates of Sec. |IV C[ these changes occur at magnetic 
length Ib ^ a on the order of the lattice scale, where 
the renormalizations are still moderate in magnitude, 
^e-ei}B / O'^w) ^ 1, i.c., for unrealistically high magnetic 
fields. For all laboratory fields, Ib/cl ^ 10 — 100 and the 
renormalizations are strong, Fc-c{Ib / o.,w) S> 1. In this 
regime, the signs and relative value of the anisotropy en- 
ergies u_L,z(i3_L) [Eq. (4.26)] cannot anymore be changed 



and, hence, the transitions cannot be induced, by varying 
the magnetic field. 

A potentially more fruitful way to induce phase tran- 
sitions could be by tilting the mag netic fiel(p2l2l!. The 
anisotropy energies uj_,z(i?j_) depend on the magnetic 
field component -B_l, perpendicular to the sample, while 
the Zeeman energy £z{B) = (J^bB depends on the to- 
tal value B = \J B\ + By . Therefore, tilting the mag- 
netic field increases the Zeeman energy tz relative to the 
anisotropy energies which makes the F phase more 
favorable. Provided at perpendicular field orientation the 
system is in a different phase, eventually, upon increasing 
B/i?_i_, the system must end up in the F phase. 

Figure [19] shows the evolution of phases as one ap- 
plies the parallel component B\\ , while keeping B±^ fixed, 
thus changing the Zeeman energy from ezi_ = fJ-sB^ to 
ez = I^-bB > ez± and keeping the anisotropy energies 
{u±,Uz) constant. It is convenient to present the phase 
diagram in the units ^f^(uj_,Uz): this way, the phase 
boundaries remain fixed, while the phase point with con- 
stant {u±, Uz) moves along the straight line from its posi- 
tion at By = to the origin at B/B± — > oo. If the system 
is in either the CAF or CDW phases at By = 0, points 
a or c, respectively, the transition occurs directly into 
the F phase upon increasing ez/ez±i via a continuous 
second-order transition in the former case and a first- 



order transition in the latter case. Remarkably, however, 
if the system is in the KD phase at -By = 0, point b, the 
transition to the F phase cannot occur directly: upon 
increasing ez, the system first makes a first-order phase 
transition into the CAF phase, and then continuously 
transitions to the F phase. 

We note, however, that, since the renormalizations of 
the isospin anisotropies are strong, such transitions may 
be quite challenging to realize. According to Sec. |IV C[ 



ahhough the bare energies uf'.^ ~ Bj_[T]K [Eq. (|4.28l) 
are comparable with the Zeeman energy ez± = }j-bJ^± 
for perpendicular field orientation, the renormalized en- 
ergies M_L_2 [Eqs. (4.29) and (4.30)] can easily exceed 
ez_L by one order of magnitude^ Therefore, achieving 
the transition by applying the parallel component By , 
would generally require large ratios B/B± > 10, i.e., 
large tilt angles. Since one should, at the same, main- 
tain a large enough perpendicular component B± to pre- 
serve the correlated quantum Hall state above the disor- 
der level {Bj_ « 3 - 5T, according to Refs. I15|17|) . the 
practical maximum of the ratio B/B± is limited by the 
maximum achievable magnetic field and disorder in the 
system. An exception concerns the special cases when, 
for perpendicular field orientation, the system is "anoma- 
lously" close to one of the KD-CAF, CDW-F, or CAF- 
F transitions lines in Fig. 16 [i.e., the point {uj_,Uz) is 



ez± away from the linejand is on the KD, CDW, 
or CAF side, respectively. In this case, smaller ratios 
B I B± ~ 1 would be sufficient for the phase transition to 
occur. 



VI. CONCLUSIONS 

In conclusion, in this paper, we studied the v = quan- 
tum Hall state in monolayer graphene in the framework 
of quantum Hall ferromagnetism, with the key empha- 
sis on the isospin anisotropies that arise from the valley 
and sublattice asymmetric short-range electron-electron 
and electron-phonon interactions. The phase diagram in 
Fig. [18] obtained in the presence of the generic isospin 
anisotropy and Zeeman effect (neglecting thermal fluc- 
tuations), consists of four phases characterized by the 
following orders: spin ferromagnetic (F), charge density 
wave (CDW), Kekule distortion (KD), and canted anti- 
ferromagnetic (CAF). To the best of our knowledge, the 
CAF phase has not been predicted before in the con- 
text of the correlated quantum Hall states in graphene. 
We took into account the Landau level mixing effects 



(Sec. IV I and found that they result in (i) the suppression 
of the stiffness ps due to screening and (ii) critical renor- 
malizations of the anisotropy energies u±^z- The latter 
has crucial implications for the physics of the i/ = state. 
First, the anisotropies are greatly enhanced and can sig- 
nificantly exceed the Zeeman energy. Second, and most 
importantly, we conclude that the short-range electron- 
electron interactions could favor any state on the generic 
phase diagram: one cannot theoretically rule out any pos- 
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sibility based just on the repulsive nature of the underly- 
ing Coulomb interactions. The leading electron-phonon 
interactions, on the other hand, always favor the Kekule 
distortion phase. 

The main open practical question is then which of the 
phases in Fig. [18] corresponds to the strongly in sulating 
i/ = state observed in transport experiment j ^^^^ ^^ l ^^ l 
For high quality suspended graphene sample^^^^, the 
two-terminal resistance was instrument-limited at i? ~ 
109 _ IQioOhm at B ~ lOT and T < IK. In the absence 
of an accurate estimate for the anisotropy energies, one 
can try to infer about the nature of the real i/ — state 
by addressing the transport properties of the phases in 
Fig.lll 

While the charged excitations o f the v = QHFM are 
gapped in the bulk for any ordeif^UsHHIlI]^ the phases 



in Fig. 18 are expected to have m arkedl y different edge 
transport behavior. Existing studieP^l^ suggest that the 
F phase has gapless edge excitations. Therefore, an ideal 
sample with the F bulk order would have a two-terminal 
resistance R = l/(2e^//i) ~ 10'*Ohm, with the fa ctor 2 
du e to two edges. On the other hand, the CDW^aMl^ 
AfUMH, and KjJ^ phases were shown to have gapped 
single-particle mean-field edge excitations, which implies 
that these phase should exhibit insulating behavior. The 
edge excitations of the CAF phase have not yet been ad- 
dressed. Note, however, that the spin-polarized F phase. 



being an isospin singlet, is the only phase of the v = Q 
QHFM that does not break the sublattice symmetry. 
This observation could be used as a "rule of thumb" argu- 
ment in speculating about the edge transport properties 
of the V = Q states. Based on this, one can expect all 
three phases, CDW, KD, and CAF, to have gapped edge 
excitations and exhibit insulating behavior, which makes 
them much stronger candidates than the F phase for the 
V — Q state realized experimentally. Distinguishing fur- 
ther between the three insulating phases based on the 
transport properties requires a more detailed analysis of 
their bulk and edge charge excitations, to be presented 
elsewherJ^U 
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